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EXECUTIVE  SUMMARY 


This  report  analyzes  the  ability  of  unintended  intercept  receivers  to 

isolate  individual  Direct  Sequence  (DS)  emitters  within  a  network  of  DS 

4 

emitters.  The  report  demonstrates  that  interception  over  1/R  propagation 
paths  requires  very  large,  ground-based  antennas  in  order  to  achieve  a  usable 
intercept  Signal-to-Noise  Ratio  (SNR).  On  the  other  hand,  free  space,  1/R^ 

i 

propagation  paths  provide  intercept  receivers  with  signal  power  levels  well 
above  tne  thermal  noise  power  even  when  the  intercept  antenna  is  quite  small. 


However,  these  small  antennas  are  not  able  to  resolve  individual  emitters  at 


typical  intercept  ranges.  Antennas  large  enough  to  resolve  individual 
emitters  are  too  large  to  put  on  aircraft  as  required  to  achieve  the  necessary 
free  space  propagation.  Hence,  the  report  finds  that  the  beamwidth  of  the 
intercept  receiver  will  typically  contain  many  DS  emitters.  Therefore  the 
receiver  bandwidth  will  contain  many  overlapping  DS  signals  each  typically 
with  a  power  level  which  is  well  above  the  thermal  noise  power  in  the 
intercept  bandwidth.  This  is  the  intercept  context  analyzed  in  this  report. 


The  analyses  therefore  confirm  that  the  many  DS  signals  mask  each  other 
at  linear  receivers  which  attempt  to  find  a  spectral  region  in  which  the  power 
of  one  signal  is  well  above  the  combined  powers  of  the  other  DS  signals.  This 
follows  because  if  the  number  of  DS  emitters  in  the  antenna  beamwidth  equals 
M+l,  then  the  SNR  for  any  one  signal  is  on  the  order  of  1/M,  which  is  expected 


to  be  low.  Hence,  to  be  a  serious  threat,  the  intercept  receiver  must  be  able 


to  detect  and  process  signals  which  are  weak  compared  to  the  total  power  in 
the  intercept  bandwidth.  Such  receivers  are  called  nonlinear  receivers. 


Nonlinear  receivers  include  the  total  power  radiometer,  the  doubler,  the 


single  channel  delay  or  autocorrelation  receiver,  the  dual  channel  or  cross 


correlation  receiver,  and  a  number  of  receiver  subsystems  called  feature 


detectors.  The  most  serious  threat  to  DS  emitters  is  found  to  be  the  cross 


correlation  receiver.  This  receiver  forms  the  cross  correlation  function  of  ,r 


the  signals  received  in  its  two  channels. 
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First  assume  that  the  cross  correlation  receiver  is  contaminated  only  by 
thermal  noise  (no  narrowband  interference).  In  this  idealized  case,  the  cross 
correlator  is  expected  to  generate  an  autocorrelation  spike  for  every  DS 
emitter  whose  location  is  such  that  the  differential  propagation  delay  of  its 
signal  between  the  two  intercept  antennas  is  more  than  just  a  few 
nanoseconds.  Each  autocorrelation  spike  can  be  converted  into  a  Line  of 
Bearing  (LOB).  However,  the  detection  of  these  autocorrelation  spikes  can  be 
greatly  complicated  by  using  signal  design  steps.  One  signal  design  step  is 
to  use  repeated,  short  composite  codes  to  produce  DS  signals  with  multiple 
spikes  in  their  autocorrelation  functions.  Common  code  segments  can  be 
inserted  into  all  DS  signals  to  produce  many  spikes  in  the  cross  correlation 
function  which  do  not  correspond  to  autocorrelation  spikes.  While  such  signal 
design  does  not  avoid  the  largest  autocorrelation  spikes,  it  does  dramatically 
increases  the  signal  processing  tasks  at  the  intercept  receiver. 

The  report  shows  that  the  major  factor  impacting  the  ability  of  the  cross 
correlator  to  detect  individual  DS  emitters  is  the  presence  of  a  dense 
environment  of  unwanted,  dynamically  varying,  narrowband  signals.  Such  is  the 
signal  environment  expected  at  airborne  interceptors.  The  report  concludes 
that  while  not  sufficient  by  themselves,  signal  design  steps  plus  operation  in 
a  dense  signal  environment  will  seriously  degrade  the  ability  of  cross 
correlators  to  resolve  and  identify  individual  DS  emitters  in  a  DS  network. 

Furthermore,  to  be  a  major  threat,  the  LOB's  from  a  DS  emitter  generated 
at  several  intercept  sites  must  be  coordinated  to  form  estimates  of  emitter 
location.  This  coordination  requires  that  the  autocorrelation  spikes  be 
"tagged"  with  a  signal  feature  as  a  means  of  identification.  Typical  tags 
which  can  be  extracted  in  the  low  SNR  environment  by  nonlinear  feature 
detectors  include  the  DS  chip  rate  and  recovery  of  the  DS  carrier.  The 
process  by  which  an  unintended  receiver  can  extract  useful  identifying  tags 
can  be  very  seriously  degraded  by  the  design  of  the  DS  signal.  For  instance 


design  procedures  which  limit  the  performance  of  chip  rate  detectors  include 
the  following: 

•  Frequency  Hopping  to  limit  the  bandwidth  of  the  bandpass  filter  to  a 
value  equal  to  about  twice  the  hopping  rate; 

•  Time  Hopping  to  limit  the  bandpass  filter  bandwidth  to  a  value  equal 
to  about  the  reciprocal  of  half  the  burst  duration; 

•  Tailoring  the  chip  shape  to  approach  a  bandlimited  spectrum  such  that 
the  spectral  line  at  the  chipping  rate  is  small  as  it  is  at  all  the 
harmonics  of  the  chipping  rate;  and 

•  Jittering  the  chip  rate  such  that  again  no  strong  spectral  line  exists 
at  the  chipping  rate. 

These  design  steps,  in  addition  to  the  fundamental  goal  of  operation  in  a 
dense  signal  environment,  should  make  the  combined  processes  of  spike 
detection  and  tag  extraction  extremely  unreliable  and  time  consuming.  In  the 
case  of  the  carrier  recovery  detector,  signal  design  steps  include  the 
following: 

•  Using  MSK  or  QPSK  spreading  to  lower  the  output  SNR  for  a  given  input 
SNR  and  TW  product; 

•  Using  digitized  band-limited  Gaussian  noise  as  spreading  sequences 
with  essentially  uniform  phase  distribution  and  Rayleigh  distributed 
amplitude; 

•  Jittering  the  carrier  frequency;  and 

•  Time  Hopping  and  Frequency  Hopping  to  limit  integration  time  as 
against  the  chip  rate  detector. 

Again  these  steps  should  be  combined  with  operation  in  a  dense  signal 
environment. 

It  is  concluded  that  using  some  of  the  above  design  steps  plus  the 
inherent  sensitivity  of  nonlinear  detectors  to  dense  environments  can  keep 
nonlinear  receivers  from  performing  adequately  even  when  excision  is  used  to 
prewhiten  the  input  and  remove  distortion  from  the  filter  output. 
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CHAPTER  1 
INTRODUCTION 

A  primary  consideration  in  denying  an  enemy  knowledge  of  friendly  force 
position  and  movement,  location  of  critical  command  centers,  and  Electronic 
Support  Measures  (ESM)  is  using  communications  signals  with  a  low  probability 
of  detection  and  interception.  One  means  of  providing  tactical  Army 
communications  with  a  Low  Probability  of  Detection  (LPD)  and  a  Low  Probability 
of  Identification  (LPI)  is  the  use  of  Direct  Sequence  (DS)  spread  spectrum 
modulation.  Not  only  does  the  wideband  DS  modulation  lower  the  signal  power 
spectral  density,  making  the  signal  difficult  to  detect  and  identify,  but  the 
signals  from  the  many  emitters  in  a  network  also  tend  to  overlap  and  mask  the 
presence  of  each  other.  The  purpose  of  the  research  reported  here  is  to 
quantitatively  assess  and  parameterize  the  LPD  and  LPI  protection  achieved  by 
a  single  DS  emitter  within  a  network  of  DS  emitters. 

The  research  reported  here  focuses  on  two  critical  emitter  detection  and 
isolation  issues.  The  first  area  investigates  the  feasibility  of  an  intercept 
receiver  capable  of  isolating  an  individual  communications  emitter  using  a 
highly  directional  antenna.  The  second  area  investigates  the  ability  of  an 
intercept  receiver  to  somewhat  limit  the  number  of  emitters  in  its  antenna 
mainlobe  and  then  to  isolate  on  an  individual  emitter  using  signal  processing 
techniques.  This  second  area  immediately  results  in  consideration  of  signal 
processing  techniques  capable  of  performing  well  when  the  power  of  the  Signal 
of  Interest  (SOI)  is  lower  than  the  sum  of  the  powers  of  the  other  signals 
present.  This  is  the  case  because  it  would  not  be  unlikely  for  the 
Signal-to-Noise  Ratio  (SNR)  for  one  emitter  among  M  emitters  in  the  antenna 
mainlobe  to  be  approximately  1/(M-1),  where  M  is  the  total  number  of  emitters 
in  the  mainlobe. 


BACKGROUND 


2.1  INTERFERENCE  IN  COMMUNICATIONS  NETWORKS 


Earlier  analysis  of  the  design  of  an  LPI  network  has  confirmed  that  it  is 
not  unusual  for  the  "noise  floor"  of  a  receiver  in  a  network  of  DS  radios  to 
be  established  by  self-interference  caused  by  the  reception  of  unwanted  DS 
signals  transmitted  by  other  emitters  in  the  network.  This  is  particularly 
the  case  in  the  low  UHF  band  (400  MHz)  where  the  propagation  attenuation  is 
lower  and  the  antennas  are  less  directional  than  at  higher  UHF  frequencies 
such  as  1800  MHz.  Mathematically  this  means  that  for  communications  power 
levels  and  ranges  within  a  ground-to-ground  network  of  direct  sequence  spread 
spectrum  radios  it  is  not  unusual  for  the  following  condition  to  prevail 
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where  the  left-hand  side  of  Equation  (2.1)  is  the  interference  power  received 

at  one  DS  radio  as  a  result  of  the  signal  power  transmitted  by  the  other 

radios  In  the  DS  network.  In  Equation  (2.1),  PT  is  the  assumed  equal 

transmitted  power  level  of  each  emitter  in  the  network,  G_  is  the  gain  of 

th  TmR 

the  antenna  of  the  m  emitter  in  the  direction  of  a  given  receiver,  G 

th  RTm 

is  the  gain  of  the  receive  antenna  in  the  direction  of  the  m  transmitter, 

/3(Rm)  is  the  propagation  attenuation  assumed  to  follow  a  1/R  dependence 

over  the  path  from  the  mfc^  transmitter  to  the  given  receiver,  W  is  the 
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spread  spectrum  bandwidth,  W  is  the  communications  baseband  bandwidth,  and 

DO 

Nq  is  the  thermal  noise  power  spectral  density  at  the  receiver.  Of 
particular  concern  is  the  attenuation  of  the  signal  from  the  nearest 
transmitter  in  the  network  whose  message  is  unwanted  at  a  given  receiver.  As 
long  as  the  conditions  indicated  by  Equation  (2.1)  prevail,  the  SNR  at  any 
receiver  is  not  significantly  increased  or  decreased  by  increases  or  decreases 
in  transmitted  power,  assuming  each  emitter  transmits  at  the  same  power.  Of 
course,  the  condition  will  not  be  satisfied  as  the  transmitted  power  is 
lowered  substantially. 
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Since  the  transmitted  power,  P^,  equals  the  transmitted  energy  per  bit. 
Eg,  times  the  bit  rate,  R^  Equation  (2.1)  may  be  rewritten  yielding  the 
following  equation: 
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If  the  interference  level  at  the  communications  receivers  could  be  lowered, 
then  the  communications  transmitted  power  could  be  lower  to  provide  more 
protection  against  unwanted  detection  and  interception.  As  noted  earlier,  one 
way  to  lower  the  interference  is  to  lower  the  transmitted  power.  However,  it 
is  not  always  possible  to  lower  the  power  enough  that  it  is  not  a  major  source 
of  interference  at  the  nearest  receiver  and  at  the  same  time  provide  reliable 
communications  at  distant  receivers  to  which  it  is  transmitting.  However,  if 
the  data  rate  is  reduced  then  the  power  could  be  reduced  by  the  same  fraction 
without  lowering  the  received  energy  per  bit.  If  the  data  rate  remains  the 
same,  then  the  power  control  issue  is  range  dependent  and  needs  to  be 
addressed  in  the  network  routing  control  scheme. 

Lowering  the  interference  could  also  be  accomplished,  and  is  accomplished 

to  a  significant  degree,  by  increasing  the  spread  spectrum  bandwidth  W  . 

bb 

With  power  control  so  range  dependent.  Equation  (2.2)  suggests  that  the 
maximum  possible  spreading  bandwidth  be  used.  The  interference  can  also  be 
lowered  by  decreasing  the  values  of  the  interference  terms  added  and  by 
reducing  the  number  of  interference  terms  added.  Lowering  the  value  of  the 
summation  may  be  achieved  by  lowering  the  number  of  emitters  in  the  network. 
This  lowers  the  number  of  terms  added,  but  more  importantly,  it  tends  to 
increase  the  range  to  the  nearest  emitter  in  the  network  thus  dramatically 
increasing  the  propagation  attenuation  /9(R)  experienced  by  the  interference. 
Hence  the  values  of  the  interference  terms  in  Equation  (2.2)  will  typically  be 
significantly  decreased  by  reducing  the  number  of  transmitters  in  the  network. 


From  the  above  discussion,  it  follows  that  both  reducing  the  data  rate  and 
the  number  of  emitters  in  the  network  will  permit  the  transmitted  power  to  be 


lowered.  However,  as  determined  below,  the  impacts  of  these  two  options  on 
the  probabilities  of  detection  and  interception  at  an  unintended  receiver  may 
be  very  different. 


2.2  INTERFERENCE  AT  UNINTENDED  RECEIVERS 


The  total  amount  of  power  lying  within  the  detection  bandwidth,  WQ,  of 
an  unintended  intercept  receiver  is  given  approximately  by 
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where  N  „W_  is  the  thermal  noise  power,  P„,  is  the  received  power  of  the 
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k  source  of  narrowband  interference, 0(R  )  is  the  propagation  attenuation 
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over  the  intercept  range  R  between  the  intercept  receiver  and  the  m  DS 
communications  emitter  and  is  the  gain  of  the  intercept  receiver  in  the 

direction  of  the  mth  emitter.  If  the  DS  emitters  and  sources  of 
interference  can  all  be  separated  or  isolated  from  one  another  using  a 
directional  antenna,  the  power  received  at  an  intercept  receiver  from  the 
mth  emitter  would  be  given  by 
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where  GJT  is  the  mainlobe  gain  of  the  highly  directional  intercept  antenna 
when  it  is  pointing  directly  at  an  individual  emitter.  In  the  situation 
suggested  by  Equation  (2.4),  it  very  likely  will  be  possible  to  detect  and 
extract  the  features  of  the  DS  signal  from  the  mth  emitter  needed  to 
intercept  or  recognize  it  as  a  signal  of  interest.  This  possibility,  and  ways 
to  reduce  the  probability  of  detection  and  intercept  when  this  case  arises, 
has  been  investigated  earlier.  The  result  was  that  increasing  power  would  not 
result  in  significantly  improved  communications,  assuming  all  emitters 
transmitted  the  same  power  level.  Rather  data  rate  or  communications  range 
had  to  be  decreased  assuming  that  the  spread  spectrum  bandwidth  had  been 
maximized  during  design  and  was  fixed  by  implementation. 


•  Ji  ">  W  'A.  V  A  A.  ^  *  A*  OWA  *.  A  A 


The  more  general  case  suggested  by  Equation  (2.3)  is  investigated  here  to 
see  how  feasible  it  is  to  field  an  intercept  receiver  capable  of  the  isolation 
needed  to  achieve  the  near  worst-case  communication  situation  depicted  by 
Equation  (2.4). 

One  option,  identified  above,  was  to  lower  the  number  of  users.  However, 
this  makes  it  more  likely  that  the  individual  emitters  remaining  in  the 
network  will  be  separated  enough  in  angle  to  permit  the  type  of  isolation 
suggested  by  Equation  (2.4).  Hence  the  advantage  of  self-masking  may  be  lost 
when  the  level  of  self-interference  is  lowered  by  reducing  the  number  of 
transmitters  in  the  network.  While  the  spread  spectrum  bandwidth,  W  ,  is 
selected  during  the  design  stage  to  be  large,  it  is  not  considered  to  be  a 
candidate  for  adaptive  variation  during  operation.  This  strongly  suggests 
that  adaptively  lowering  the  data  rate  may  be  a  much  preferred  way  of  lowering 
transmitted  power  to  achieve  LPD  and  LPI.  This  lowers  the  SNR  at  the 
intercept  reciever  and  also  has  the  potential  of  retaining  the  very  desirable 
mutually  self-masking  aspects  of  the  emitters  within  a  network  of  DS  radios. 
This  possibility,  and  factors  related  to  it,  will  be  carefully  studied  in  the 
next  section  of  this  report. 


CHAPTER  3 
TECHNICAL  RESULTS 


3 . 1  INTRODUCTION 


This  section  of  the  report  is  divided  into  two  major  subsections. 
Subsection  3.2  analyzes  the  antenna  issues  relevant  to  the  detection  and 
isolation  of  individual  emitters.  This  subsection  derives  the  fact  that  at 
UHF  frequencies  near  400  MHz  that  antennas  capable  of  isolating  individual 
emitters  in  typical  tactical  geometries  are  prohibitively  large.  At  1800  MHz 
antennas  capable  of  emitter  isolation  are  unlikely,  but  possible.  Subsection 
3.3  considers  the  case  when  the  intercept  antenna  is  able  to  reduce  the  number 
of  emitters  in  its  mainlobe  to  a  number  M+l,  where  M  is  significantly  larger 
than  unity.  Further  isolation  down  to  individual  emitters  would  then  depend 
on  extracting  an  emitter  feature  or  features  by  which  the  individual  emitters 
could  be  identified. 

3.2  ANTENNA  ISSUES  IN  EMITTER  DETECTION  AND  ISOLATION 


The  factors  which  most  influence  the  ability  of  an  unintended  receiver  to 
detect  and  isolate  an  individual  emitter  within  a  network  of  DS  emitters  are 
as  follows: 

a.  the  relative  communications  and  intercept  ranges; 

b.  the  propagation  over  the  communications  and  intercept  paths; 

c.  the  gain  and  beamwidth  of  the  intercept  antenna;  and 

d.  the  intercept  bandwidth  and  signal  processing  strategy. 

The  first  three  of  these  issues  will  be  addressed  in  this  subsection  and 
the  fourth  issue  will  be  the  focus  of  the  next  subsection.  The  approach  taken 
in  this  subsection  is  to  determine  the  relationships  between  the  first  three 
factors  which  result  in  the  ability  to  detect  and  isolate  individual  emitters 
within  a  network  of  DS  emitters.  The  realizability  or  ability  to  actually 
field  an  intercept  system  which  embodies  these  factors  will  then  be  assessed. 
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The  goal  will  be  to  see  if  any  realizable  combination  of  the  above  factors  can 
be  found  which  would  permit  individual  emitters  to  be  detected  and  recognized. 

The  three  factors  addressed  in  this  subsection  will  be  analyzed  as 
follows.  The  intercept  range  expected  will  lead  to  an  indication  of  the 
antenna  gain  necessary  to  achieve  an  adequate  intercept  Signal-to-Noise  Ratio 
(SNR).  Focus  will  also  be  placed  on  the  minimum  angular  separation  expected 
between  emitters  (as  seen  by  the  interceptor)  for  typical  communications  and 
intercept  ranges.  This  angle  will  provide  an  indication  of  the  antenna 
beamwidth  needed  to  resolve  or  isolate  the  emitters  in  angle.  Assuming 
communications  center  operating  frequencies  of  400  and  1800  MHz,  both  the  gain 
required  for  adequate  intercept  SNR  and  the  beamwidth  required  for  emitter 
isolation  lead  directly  to  required  values  of  the  diameter  of  a  parabolic 
intercept  antenna  dish.  The  size  of  the  antenna  dish  required  will  be  used  as 
the  principal  measure  of  the  realizability  or  fieldability  of  a  candidate 
intercept  system.  The  comparison  between  the  antenna  size  needed  to  isolate 
emitters  and  the  size  needed  to  achieve  a  required  intercept  SNR  is  critical 
in  determining  the  relative  risk  of  being  isolated  in  angle  or  being  detected 
in  a  dense  signal  environment. 

The  second  important  factor  is  propagation.  The  propagation  path  for  the 

4 

communications  will  be  assumed  to  be  1/R  in  all  cases.  The  intercept 
propagation  path  will  typically  be  assumed  to  be  1/R^,  but  the  impact  of  a 

4 

1/R  intercept  path  will  be  noted. 

Two.  major  detection  strategies  will  be  investigated.  The  first  will  be 
linear  Fourier  power  spectral  analysis  as  a  means  of  detecting  the  DS  signal 
in  the  noise  and  interference  environment.  The  ability  of  the  antenna  to 
supply  an  adequate  SNR  will  be  a  key  issue.  The  second  detection  strategy,  to 
be  addressed  in  Subsection  3.3,  will  be  the  use  of  a  nonlinear  cross 
correlation  receiver.  In  this  case,  the  ability  of  the  antenna  to  limit  the 
number  of  other  emitters  and  sources  of  nonwhite  interference  power  in  the 
mainlobe  will  be  critical. 


3-2 


The  first  issue  to  be  investigated  is  the  size  of  antennas  needed  to 
isolate  or  resolve  individual  DS  emitters  within  a  network  of  DS  emitters. 

This  analysis  will  assume  that  the  unintended  receiver  uses  a  parabolic  dish 
or  that  the  antenna  reflector  is  at  least  parabolic  in  the  azimuth  plane. 

Under  these  assumptions,  for  a  parabolic  dish  with  an  illumination  factor  of 
0.54,  the  gain  relative  to  an  isotropic  antenna  in  dB  is  given  by  [1] 

CdB  =  20  log  fmz  *  20  log  Dft  -  52.6  (3.1) 

where  f  is  the  frequency  in  megahertz  and  D_.  is  the  dish  diameter  (or 
MHZ  £t 

aperture  in  the  azimuth  plane)  in  feet. 

Two  frequencies  are  of  particular  interest.  These  are  400  MHz  and  1800 
MHz.  At  these  frequencies  the  gain  in  dB,  from  Equation  (1)  becomes 


GdB  =  20  log  Dffc  -  0.56  at  400  MHz  (3.2a 

=  20  log  Dft  +  12.51  at  1800  MHz  (3.2t 

The  ability  of  an  antenna  to  resolve  an  emitter  depends  on  its  3  dB 
beamwidth.  Generally,  in  a  given  geometry,  emitters  lying  within  the  3  dB 
beamwidth  of  an  antenna  are  not  resolvable  in  the  given  geometric  situation 
using  that  antenna.  The  3  dB  beamwidth,  in  radians,  of  a  parabolic  dish 
antenna  is  given  by  [1] 

(3.3) 


at  400  MHz  (3.4a 

at  1800  MHz  (3.4b 


[1]  Reference  Data  for  Radio  Engineers,  ITT  Corp.,  Fourth  Edition,  1949,  p.753 
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1.223  x  10' 
fMHz°ft 


which  yields 

a  3.06 


e 


“ft 

0.679 


ft 


Using  Equations  (3.2)  and  (3.4),  the  relationship  between  the  gain  in  dB  and 
the  3  dB  beamwidth  (resolution)  in  radians  is  given  by 


0dB  .  9.14  -  20  log 


(3.5 


Two  approaches  will  be  used  to  predict  how  small  the  antenna  resolution 

must  be  in  order  to  isolate  individual  emitters.  The  first  approach,  as  used 

in  the  initial  DS  network  investigation  at  E-Systems,  assumes  a  network  of  M 

emitters  is  active  over  an  area  A.  The  average  communications  distance  (d  ) 

o 

is  then  defined  as  the  diameter  of  the  circle  whose  area  is  such  that 


A 


From  Equation  (6),  the  average  communication  distance  is  given  by 


(3.6 


(3.7 


It  follows  that  the  antenna  3  dB  beamwidth  should  satisfy  the  relationship 

<o 


Ae~ 


(3.8 


where  Rj  is  the  intercept  range.  To  insure  that  consistent  units  are  used. 


A  will  be  given  in  square  miles. 


in  miles,  and  Rj  will  be  in  miles. 


Using  Equation  (3.4)  in  (3.8)  yields 


Ri  = 

where  Cf  equals  3.06  at  400  MHz  and  0.679  at  1800  MHz. 

Equation  (3.9)  is  used  to  plot  the  antenna  diameter,  Dffc,  as  a  function 
of  the  intercept  range,  R  ,  for  the  case  when  A  equals  100  square  miles  and 
M  equals  10,  50  and  100.  The  result  when  the  communication  center  frequency 
equals  400  MHz  is  shown  in  Figure  3-1.  For  an  1800  MHz  center  frequency,  the 
results  are  shown  in  Figure  3-2.  From  these  figures,  it  is  clear  that 
emitters  at  1800  MHz  are  much  more  easily  resolved  than  emitters  at  400  MHz. 
Based  on  the  "average-communications-range"  model,  it  appears  that  50  to  100 
users  at  400  MHz  distributed  over  100  square  miles  would  be  resolved  by  an 


2  A 


Vm 


(3. 


(miles) 


intercepter  at  a  range  of  more  than  10  miles  only  if  the  antenna  diameter  were 
greater  than  about  25  feet.  The  same  25  foot  antenna  would  provide  isolation 
out  to  about  50  miles  at  1800  MHz. 


The  fact  that  a  25  foot  antenna  would  almost  surely  not  be  airborne  will 

4  2 

be  considered  later  when  the  impact  of  1/R  versus  free  space  1/R 
propagation  is  analyzed. 


The  second,  more  direct,  approach  to  intercept  antenna  requirements  is  to 
set  the  antenna  3  dB  beamwidth  equal  to  the  separation  between  emitters,  R^, 
divided  by  the  intercept  range,  R  .  In  this  case.  Equation  (3.4)  yields 


**• 


(3.10) 


Equation  (3.10)  is  used  to  plot  as  a  function  of  R^  for  various 

values  of  communication  range,  R  .  The  results  for  400  MHz  communications 

c 

is  shown  in  Figure  3-3  and  for  1800  MHz  in  Figure  3-4.  From  Figure  3-3,  it  is 
apparent  that  very  sizable  antennas  are  required  for  satisfactory  resolution 
at  400  MHz.  Note  that  if  the  communications  link  is  not  perpendicular  to  the 
line  of  intercept,  even  larger  antennas  are  required.  Hence,  the  antenna  size 
plotted  is  the  smallest  capable  of  isolating  the  emitters  with  the  specified 
separation.  Figure  3-4  shows  that  an  intercept  range  of  less  than  about  20 
miles  is  required  if  emitters  at  1800  MHz  separated  by  one-half  mile  are  to  be 
resolved  with  a  25  foot  antenna. 

Summarizing,  Figures  3-1  through  3-4  reveal  that  to  isolate  communications 
centered  around  400  MHz  requires  either  very  large  antennas  (greater  than  25 
foot  diameter)  or  very  close  intercept  locations  (less  than  5  or  10  miles). 
Antennas  this  large  would  not  be  mobile  but  could  exist  as  fixed-site 
installations.  For  fixed  sites,  the  intercept  range  is  expected  to  increase 
resulting  in  even  larger  antennas.  Hence,  isolation  of  emitters  at  400  MHz 
seems  infeasible  for  most  tactical  communications  emitters  at  1800  MHz  appears 
feasible,  but  not  likely,  except  at  intercept  ranges  less  than  about  10-20 
miles  for  emitters  separated  by  one-half  mile.  It  appears  that  a  10  foot 


r 
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S  N 


I , miles 


antenna,  consistent  with  airborne  intercept,  would  require  an  intercept  range 
of  7,5  miles  to  resolve  emitters  separated  by  one-half  mile  and  15  miles  for 
emitters  separated  by  one  mile.  It  should  be  noted  that  propagation  and 
signal  strength  considerations  as  discussed  in  the  next  section  relate 
directly  to  the  ability  of  a  receiver  to  isolate  DS  emitters. 

The  second  issue  critical  to  emitter  isolation  and  detection  is  the 
intercept  Signal-to-Noise  Ratio  (SNR).  The  SNR,  in  turn,  is  strongly 
dependent  on  the  propagation  mode  between  the  emitters  and  the  intercept 
receiver.  These  issues  are  discussed  in  this  section  of  the  report. 


A  DS  emitter  spread  its  transmitted  power,  P^,  over  its  spread  spectrum 


bandwidth,  W 


ss 


Assuming  that  the  transmitted  signal  power  spectral  equals 


P  /W  ,  the  power  in  an  intercept  bandwidth,  W  ,  is  given  by 
X  S  S  1 


Q  = 


P  (W  /W  )  G„ 

T  I  ss  TI  IT 


ai(RI)LI 


(3.1 


where  is  the  gain  of  the  communications  transmit  antenna  in  the 


direction  of  the  unintended  intercept  receiver,  GIT  is  the  gain  of  the 


intercept  antenna  in  the  direction  of  the  communications  emitter,  a^R^) 


is  the  propagation  attenuation  over  the  intercept  path  of  range  Rj,  and  L ^ 
is  the  loss  in  the  intercept  receiver.  Assuming  the  intercept  receiver  is 


thermal  noise  limited,  at  noise  power  spectral  density  NoI»  the  intercept 


SNR,  denoted  by  Q/Nj  is  given  by 


Q_ 

N 


P  G  G 
T  HIT 


aT(RT)NrtTW  L. 
I  I  01  ss  I 


(3.1 


since  the  intercept  thermal  noise  is  equal  to  N  W^.  In  order  for  Q/Nj 
to  exceed  the  required  value,  denoted  by  (Q/Nj)*,  Equation  (3.12)  shows 
the  intercept  antenna  gain  must  satisfy  the  relationship 


However,  the  gain  of  a  parabolic  antenna  can  be  written  in  the  form 


from  which  Equation  (3.1)  follows.  Substituting  Equation  (14)  into  Equation 
(3.13)  yields 


Dft>2.87  x  10 


V 


(Q/N_)*W  N  L 
I  ss  ol  I 

p""g  [G  /G  T 
T  II  1  TI  TRJ 


I, mi 


(3.15 


when  the  intercept  path  is  free-space  and  when  Rj  is  the  intercept  range 

in  miles.  Equation  (3.15)  is  used  to  plot  D  versus  R_  .  from  typical 

ft  I, mi 

values  of  the  parameters  under  the  radical  sign.  The  assumed  values  are 

(Q/Nj)*  equals  10,  P^G,^  the  transmitted  Effective  Radiated  Power  (ERP) 

is  10  watts,  N  _  is  -196  dBW/Hz,  L_  is  3  dB  and  the  gain  ratio 
ol  I 

(G  /G  ]  equals  unity.  The  antenna  diameter  required  to  achieve  the 

X  X  TK 

necessary  value  of  (Q/Nj)*  for  three  values  of  spread  spectrum  bandwidth  are 
shown  in  Figure  3-5.  Note  that  the  diameter  is  quite  small,  relative  to  the 
diameter  needed  to  resolve  emitters,  even  with  the  reasonably  low  transmitted 
ERP  of  10  watts.  From  Figure  3-5,  it  is  clear  that  for  free  space  propagation 
the  size  of  the  antenna  required  to  achieve  useful  values  of  intercept  SNR  is 
small  compared  to  the  size  needed  to  resolve  emitters.  The  size  is  small  even 
at  ranges  compatible  with  airborne  intercept  from  many  tens  of  miles. 


when  the  intercept  propagation  path  is  not  free  space,  the  attenuation  can 
be  written  in  the  form 

a(RT)  =  C ( f , R  )a  ( f , R  ) (R  /R  )4  (3.16 

1  O  £3  O  I  O 

4 

for  typical  1/R  ,  lower  UHF  ground-to-ground  communications  over  intercept 

ranges  of  about  0.1  miles  to  10  miles.  In  Equation  (3.16)  ,  C(f,RQ)  is  a 

constant  which,  depending  on  the  heights  of  antennas,  could  reasonably  equal 

35  dB  at  a  reference  range,  R  ,  of  1  mile.  The  factor  a.  (f,R  )  is 

2  o  f  s  o 

the  free-space,  1/R  attenuation  at  the  communications  center  frequency  and 
reference  range  of  1  mile.  Therefore,  it  follows  that. 


Figure  3-5.  Antenna  Diameter  Required  to  Achieve  Useful  Intercept  SNR  Against 
Signal  with  Indicates  Spread  Spectrum  Bandwidth  (Free  Space 
Intercept  Assumed). 


=  [36.58  20  log  f 


(3.17) 


a£S<f'V  * 


L-^J 


MHz 


+  20 


log  R  . ]  dB 
o ,  mi 


-  C36-58  *  20  l0*  WdB 

at  the  referene  range  of  1  mile. 


In,  dB, 


the  1/R4 


attenuation  is  therefore 


a(R  )  =  71.58  +  20  log  f  ♦  40  log  R, 

I  ’  MHz  ’  I, mi 


(3.18) 


for  R  .  between  about  0.1  and  10  miles.  Using  Equation  (3.17),  the 
I  $  mi 

diameter  of  the  intercept  antenna  becomes 


Dft>2. 
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(3.19) 


Using  C(f,l)  equal  to  35  dB,  Equation  (3.19)  is  used  to  plot  versus  the 

intercept  range.  Note  the  dependence  on  the  square  of  the  intercept  range. 

For  the  same  values  of  the  parameters  under  the  radical  sign.  Figure  3-6  is  a 

plot  of  versus  R^  for  three  values  of  the  spread  spectrum  bandwidth. 

From  Figure  3-6,  note  that  the  antenna  diameter  required  is  much  larger  for 
4  2 

1/R  than  1/R  propagation.  For  instance  an  antenna  diameter  of  more  than 
25  feet  is  required  to  intercept  a  signal  with  a  10  MHz  spread  bandwidth  from 
an  intercept  range  of  only  5  miles  under  the  assumed  conditions. 


2 

From  Figures  3-5  and  3-6,  the  conclusion  is  that  for  free-space,  1/R 
propagation  the  intercept  antenna  diameter  is  quite  small  and  easily 
realized.  The  antenna  is  much  larger  for  ground-to-ground  intercept  of  OS, 
UHF  signals. 


Overall,  the  conclusions  from  Figures  3-1  through  3-6  are  that  while 
resolving  emitters  at  1800  MHz  is  feasible,  it  is  most  likely  that  intercept 
antennas  will  be  able  to  achieve  intercepted  signal  power  well  above  thermal 
noise  levels  but  that  several  emitters  would  typically  be  within  the  antenna 
3  dB  beamwidth,  particularly  for  communications  at  400  MHz. 


Based  on  Figure  3-6,  1/R  propagation  will  be  assumed  for  the  intercept 
path.  From  Figure  3-5,  a  25  foot  antenna  will  provide  an  intercepted  signal 
power  well  above  the  thermal  noise  floor  when  the  intercept  propagation  path 
is  1/R^.  Using  Equation  (3.10)  and  solving  for  the  communication  range  with 
a  25  foot  antenna  yields  the  following  condition  under  which  the 
communications  emitters  are  not  resolved 


Equation  (20)  yields  Figure  3-7,  which  plots  the  minimum  resolvable 
emitter  separation  using  a  25  foot  antenna  as  a  function  of  intercept  range. 
The  figure  suggests  that  many  emitters  will  be  within  the  3  dB  beamwidth  of  an 
intercept  antenna,  particularly  at  400  MHz.  For  instance,  at  an  intercept 
range  of  40  miles,  the  minimum  resolvable  separation  is  about  5  miles  at  400 
MHz.  A  circle  with  a  diameter  of  5  miles  has  an  area  of  almost  20  square 
miles.  Using  Equation  (3.7),  the  value  of  M  such  that  the  average 
communications  range,  {  dQ),  e<Jua*s  one-half  mile  is  101.  Hence  many 
emitters  are  expected  to  be  in  the  3  dB  beamwidth  and  potentially  mask  the 
presence  of  one  another. 

The  ability  to  extract  a  feature  which  could  allow  the  individual  emitters 
to  be  identified  is  discussed  in  the  next  section. 

3.3  EMITTER  DETECTION  AND  ISOLATION  USING  NON-LINEAR  RECEIVER 

3.3.1  Background  to  Non-Linear  Intercept 

It  may  be  somewhat  surprising  that  DS  signals  which  are  expected  to  have 
such  low  spectral  density  do  in  fact  arrive  at  receivers  well  above  the  noise 
power  spectral  density.  At  friendly  receivers,  to  which  the  DS  signal  is  not 
being  transmitted,  this  causes  network  interference.  At  unfriendly  intercept 
receivers,  the  DS  signal  is  well  above  the  noise  and  provides  an  appreciable 
intercept  SNR.  These  issues  will  now  be  investigated  as  preliminary  to 
determining  when  in  fact  non-linear  receivers  are  needed  because  the  signal  is 
below  the  receiver  noise. 
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Figure  3-7.  Minimum  Resolvable  (Cross  Beam)  Separation  for  a  25  foot  Antenna 
Versus  the  Intercept  Range. 
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It  would  not  be  unusual  for  a  UHF  ground-to-ground  (1/R  )  communications 


system  to  transmit  on  ERP^  of  10  to  20  dBW.  This  power,  when  spread 


uniformly  over  the  spread  spectrum  bandwidth  W  ,  yields  an  intercepted 

s  s 


power  spectral  density  given  by 
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(3.21) 
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where  1/R  intercept  propagation  is  assumed  and  the  variables  are  as  defined 
in  the  previous  section.  Assuming  the  intercept  thermal  noise  density  is 


N  .  the  ratio  Q  /N  _  in  dB  may  be  written  as  follows: 
ol  o  ol  z 


^o^ol’dB  *  ( e®pt } dBW  *  [GTR/GTI]dB  +  (GIT}dB 


(3.22) 
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Equation  (3.2)  can  be  used  to  calculate  (G  )  and  Equation  (3.17)  yields 

IT  Go 


the  free-space  attenuation.  Assuming  a  loss  of  6  dB  and  an  8  dB  noise  figure 
Equation  (3.22)  becomes 


‘WdB  *  100-S  *  <E8PT>dBW  *  2<Dft>dBft  -  ‘WdB 


"^ss^dBHz  "  2(RI,mi)dBmi 


(3.23) 


As  an  example,  suppose  the  ERP^  is  10  dBW,  the  intercept  antenna  has  an  8 


foot  diameter,  and  the  mainlobe-to-sidelobe  ratio  is  15  dB.  Then  Equation 


(3.23)  becomes 


(Q  /N  =  113.8  -  (W  )._„  -  2 (R  ,).D  , 

o  ol  dB  ss  dBHz  I, mi  dBmi 


(3.24) 
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The  signal  power  spectral  density  is  10  dB  or  more  greater  than  the  noise 
spectral  density  whenever 


(W  > ♦  2(Rt  • •< 103.8 
ss  dBHz  I, mi  dBmi  — 


(3.25) 


Even  with  a  spread  spectrum  bandwidth  of  100  MHz  (80  dBHz),  this  means  that  at 


intercept  ranges  less  than  15.5  miles  the  signal  is  well  above  the  intercept 
thermal  noise.  The  noise  floor  of  a  friendly  receiver  to  which  the 
propagation  path  is  free  space,  but  to  which  the  signal  is  not  intended,  would 
likewise  be  dominated  by  the  DS  signal  rather  than  the  thermal  noise. 


To  determine  the  impact  of  the  anticipated  1/R  propagation  path  between 
friendly  receivers.  Equation  (3.21)  is  rewritten  to  replace  the  free-space 


attenuation  with  a(R^)  given  by  Equation  (3.16).  In  this  case  Rj  would 


represent  the  range  to  the  interfering  communications  transmitter.  At  400 


MHz,  and  using  C(f,f  )  as  35  dB  at  one  mile,  we  obtain 
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Using  the  same  values  as  follow  Equation  (3.23),  Equation  (3.26)  becomes 


(Q  /N  _)._  =  78.8  -  (W  )._„  -  4  ( R  .  ).„  . 

o  ol  dB  ss  dBHz  I, mi  dBmi 


(3.27) 


It  follows  that  the  noise  floor  would  be  increased  3  dB  or  more  because  of 
the  addition  of  the  signal  power  density  whenever 


(W  )  ._  +  4 ( R  . )  .  <  78.8 

ss  dB  I, mi  dBmi 


(3.28) 


With  a  spread  spectrum  bandwidth  of  80  dBHz,  this  occurs  for  an 
interference  range  of  0.93  miles  or  less.  Hence  friendly  transmitters  closer 


than  about  a  mile  are  potential  sources  of  significant  self-interference  for 
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ground-to-ground  1/R  communications  at  400  MHz.  Many  transmitters  simply 
increase  the  potential  for  self-interference .  Interference  effects  could 


apparently  be  improved  by  increasing  frequency  to  1800  MHz  thus  achieving 
greater  attenuation,  higher  mainlobe-to-sidelobe  ratios  and  larger  spreading 
bandwidth.  However,  as  shown  earlier,  intercept  at  1800  MHz  requires  a 
significantly  smaller  intercept  antenna  to  achieve  isolation  of  communication 
emitters.  Multiple  mainlobe  emitters  may  well  be  needed  to  complicate  the 
operations  at  unintended  interceptors. 

From  Equation  (3.25),  it  is  apparent  that  the  spread  spectrum  signal  power 
in  an  appropriate  intercept  bandwidth  may  well  be  substantially  above  the 
thermal  noise  power  in  a  variety  of  intercept  situations.  This  being  the 
case,  the  remainder  of  this  section  will  assume  that  the  intercepted  power  is 
well  above  the  thermal  noise.  From  the  Low  Probability  of  Detection  (LPD) 
viewpoint,  this  is  conservative  since  whatever  contribution  the  noise  makes 
will  be  to  degrade  the  performance  of  detection  and  feature  extraction 
receivers . 

3.2.2  Detection  Using  a  Dual  Channel  Cross  Correlation  Receiver 

The  dual  channel  cross  correlation  receiver  is  one  type  of  non-linear 
receiver  which  in  theory  is  capable  of  detecting  signals  even  when  the  signal 
power  spectral  density  is  below  the  thermal  noise  power  spectral  density  at 
all  frequencies.  A  functional  block  diagram  of  a  cross  correlation  receiver 
is  shown  in  Figure  3-8.  The  signals  from  channels  A  and  B  provide  a 
multiplicative  reference  for  each  other.  When  the  signal  is  much  stronger 
than  the  noise,  the  performance  is  similar  to  a  matched-filter  receiver  with  a 
noisy,  rather  than  perfect,  local  replica  of  the  incoming  signal  to  be 
detected.  When  the  intercepted  signal  power,  Q,  is  well  below  the  intercept 
power,  N  ,  the  output  signal  to  noise  ratio,  (SNR)^,  is  given  by 

(SNR)  =  2 ( Q/N) \  TW  (3.29 

o  I 

This  expression  is  typical  of  the  output  SNR  of  squaring  or  correlation 
receivers.  In  other  receivers,  the  factor  of  two  is  replaced  by  another 
number  less  than  two  and  typically  between  1  and  1/4.  The  factor  is  two  in 
Equation  (3.29)  because  of  the  two  antennas  used  by  the  cross  correlation 


receiver.  Equation  (3.29)  suggests  that  provided  the  input  bandwidth  W 
contains  some  signal  power  Q  that  integration  over  a  sufficiently  long 
interval  T  will  result  in  a  useful  output  SNR  regardless  of  how  small  the 
input  SNR  denoted  by  (Q/Nj)  becomes. 

In  practice  the  processing  time  T  is  limited  by  a  number  of  issues 
including  receiver  mismatch,  gain  fluctuations,  signal  variations,  and 
receiver  precision  and  dynamic  range.  The  ability  of  the  cross  correlation 
receiver  to  detect  signals  below  the  interference  level  is  very  much  dependent 
on  the  strength,  amplitude  variation,  number,  and  dynamic  behavior  of  these 
nonwhite  sources  of  interference.  Correlated  sources  of  interference  entering 
both  channels  also  degrade  receiver  detection  performance. 

The  antenna  beamwidth  of  the  correlation  viewer  is  assumed  to  be  quite 
broad  so  that  special  antenna  tracking  is  not  required.  Hence  the  signals 
from  many  emitters  in  the  DS  network  are  expected  to  enter  the  receiver. 

These  DS  signals  are  expressed  in  the  from 

s  (t)  =  \/2S  b  (t)u  (t)  cos  [2irf  t  +  0  ]  (3.30) 

m  m  m  m  cm 

where  b  (t)  is  the  bit  phase  of  the  mth  emitter  at  time  t;  u  (t)  is  the 
m  ^  ^ 

chip  stream  of  the  m  emitter  (which  is  assumed  to  be  statistically 
independent  from  emitter  to  emitter);  f c  is  the  center  frequency  assumed  to 

be  essentially  the  same  for  all  emitters  in  the  network;  and  0  is  the 

th  ^ 

carrier  phase  of  the  m  carrier  relative  to  0q,  which  is  assumed  to  be 

zero.  This  means  that  the  signal  component  of  the  composite  waveform  entering 

the  multiplier  in  Figure  3.8  from  channel  A  is  given  by 

s  ft-A)  =V2S~  b  (t-A)u  ( t-A) cos  [2fff  (t-A)+0  ]  (2 

mA  m  m  m  cm 

The  signal  component  from  channel  B  is  given  by 

s  (t-x-5  )  =  \J 2S  b  (t-T-6  )u  (t-T-6  ) 
m  m  mm  mm  m 

cos  [2irf  ( t-T-6  )+  0  )  (3.32) 

c  mm 
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where  is  the  differential  propagation  delay  of  the  mth  signal  between 
the  two  antennas . 

The  threat  posed  by  the  cross  correlator  is  that  by  proper  selection  of 

the  variable  delay  x  in  channel  B  to  a  particular  value,  x  *,  the  mth 

in 

signal  in  channels  A  and  B  will  be  coherent  resulting  in  an  output  of  the 

integrator  equal  to  S  plus  distortion.  This  particular  value  of  x  * 

th  ^  ^ 

which  brings  the  m  signal  into  coincidence  at  the  multiplier  is  equal  to 

A  -  where  A  is  the  maximum  differential  delay  equal  to  at  least 

the  separation  between  antenna  phase  centers  divided  by  the  speed  of  light. 

The  most  negative  propagation  delay  occurs  when  the  signal  arrives  end-fire 

from  antenna  B  to  antenna  A.  Hence,  the  maximum  required  variable  signal 

delay  equal  2A. 

Looking  at  Equations  (31)  and  (32)  it  is  clear  that  the  mean  shape  of  the 

integrator  output  near  x^*  is  determined  by  te  autocorrelation  function  of 

u  (t).  Since  this  function  is  triangular  spike  band-limited  to  the  receiver 
m 

bandwidth  W,  the  duration  of  the  spike  is  approximately  1/W.  However,  since  w 

is  very  nearly  the  spread  spectrum  bandwidth,  W  ,  this  spike  is  very  narrow 

s  s 

having  a  duration  about  equal  to  the  chip  interval.  The  sharpness  of  the 

autocorrelation  spike  also  dictates  the  increment  by  which  the  variable  delay 

should  be  changed.  Assuming  a  delay  increment  which  is  a  small  fraction,  say 

1/5  or  1/10,  of  the  chip  interval  the  output  of  the  integrator  as  a  function 

of  the  delay  x  is  a  sequence  of  autocorrelation  spikes.  The  delays  x  * 

th  ^ 

at  which  they  occur  defines  a  Line  of  Bearing  (LOB)  on  which  the  m  emitter 
is  located.  If  an  identifying  feature  of  the  signal  transmitted  by  the  mth 
emitter  an  be  extracted  by  which  LOB's  can  be  exchanged  and  combined  among 
intercept  receivers,  then  the  m^  emitter  can  be  located.  One  such  tag  or 
identifying  feature  is  the  precise  chip  rate  of  the  mth  emitter.  Other  tags 
include  the  exact  center  frequency,  corrected  for  doppler  if  necessary, 
received  power  level  in  some  situations,  signal  timing  as  the  transmitter  is 
tuned  off  and  on,  and  perhaps  certain  patterns  that  the  autocorrelation  spikes 
produce.  Tags  of  this  type  may  pose  a  very  severe  threat  to  DS  emitters  when 
used  in  conjunction  with  LOB  information  derived  from  a  cross  correlation 


receiver.  Various  tags  will  be  discussed  later  in  this  section,  particularly 
the  chip  rate  and  the  center  or  carrier  frequency  of  a  DS  signal. 


The  critical  issue  is  if  a  tag  exists  how  can  it  be  unambiguously 
associated  with  a  certain  autocorrelation  spike  and  hence  a  certain  LOB.  It 
appears  that  once  is  carefully  determined  based  on  the  location  of  the 
autocorrelation  spike,  then  additional  feature  detectors  can  be  switched  into 
the  receiver.  For  instance  an  addition  delay  in  one  channel  of  approximately 
half-a-chip  interval  could  yield  the  input  to  a  chip  rate  detector.  The 
square  of  the  mth  signal  achieved  when  is  used  could  provide 
information  on  the  carrier  center  frequency. 

The  correlation  spike  of  height  is  distorted  by  the  effects  of  noise, 
the  other  DS  signals,  and  other  interfering  signals.  When  taking  these 
sources  into  account,  the  output  SNR  for  signal  sQ(t)  as  the  square  of  the 
expected  value  divided  by  the  variance  of  the  integrator  output,  can  be 
written  in  the  form 


(SNR) 


D  +  D__  +  D„_  +  D.,_  +  D  ♦  D 
SS  SP  NS  NP  NN  PP 


(3.3 


where  the  particular  delay  is  assumed  to  exist  in  the  cross 

correlation  receiver.  Th  term  D  is  the  distortion  power  due  to  the 

ss 

multiple  signals  interacting  with  one  another,  Dgp  is  the  distortion  power 
of  the  signals  interacting  with  narrowband  interferers,  D  is  the  result  of 
noise  interacting  with  the  signals,  DNp  is  the  power  caused  by  the 
interaction  of  noise  and  the  narrowband  interferers  and  and  Dpp  are 
the  powers  caused  by  the  interaction  of  noise  with  itself  and  narrowband 
inteference  with  itself  respectively. 


The  mathematical  expressions  for  each  of  the  distortion  power  terms  are 
given  below: 
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where  the  value  1.36  arises  from  the  matching  spectral  taper  of  the 
interfering  DS  signals  and  L  is  the  number  of  signals  which  interfere  with 


signal  s  ( t) 
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where  M  is  the  number  of  narrowband  inteferers  with  powers  P  and  the 

2  m 
[(sin  x)/x]  factor  accounts  for  the  offset  of  the  narrowband  interferers 


from  the  center  DS  frequency. 
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where  N  is  the  thermal  noise  power  spectral  density  of  the  cross 


correlation  receiver ,  assumed  the  same  in  both  channels. 
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where  N^  equals  the  intercept  thermal  noise  power  in  the  intercept  bandwidth 


N  2W 
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Of  special  interest  in  a  test  case  is  the  situation  when  P^  equals  zero 


for  all  k  and  there  is  only  one  signal,  s  (t).  In  this  case  Equations 

o 


(3.33)  to  (3.39)  yield 
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which  is  a  known  result  and  agrees  with  Equation  (3.29)  when  the  input  SNR/ 
denoted  by  (Q/N^),  is  much  smaller  than  unity. 

As  indicated  earlier,  of  particular  interest  is  the  case  when  the  signal 

power  is  well  above  the  thermal  noise  power.  In  this  case,  the  contributions 

of  D  and  D  can  be  ignored  and  the  output  SNR  becomes 
NS  NN 


(SNR)  =  - - — - - - rr - 

o  DSs  *  DSP  +  dNP  *  DPP 


(3.41) 


where  the  D  term  is  retained  because  the  narrowband  interferers  can  be  so 
NP 

strong  compared  to  the  signals. 

Two  cases  related  to  Equation  (3.41)  will  be  considered.  First,  is  the 
case  of  no  narrowband  interference  for  which 


(SNR) 
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When  all  powers  are  the  same  Equation  (3.42)  yields 


(SNR) 
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which  generally  speaking  is  expected  to  be  large.  However,  if  S^  is  small 
compared  to  some  Sm>  m  not  equal  to  zero,  there  is  an  opportunity  for  self 
masking  even  in  the  cross  correlator.  Assume,  for  instance,  that  is  much 
stronger  than  all  other  signals.  Then  Equation  (3.42)  yields 
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=  1 . 36 ( S  /S, )  TW 

o  1  ss 

which  shows  how  one  strong  signal  can  mask  a  much  weaker  one.  However,  the 

output  SNR  for  emitter  number  "1"  would  be  very  large.  This  could  result  in 

emitter  number  "1"  either  lowering  its  power  or  perhaps  being  destroyed.  This 

would  then  allow  the  TW  factor  in  the  output  SNR  to  yield  an  appreciable 

s  s 

output  SNR.  It  appears  that  strong  signals,  and  their  emitters,  could  be 
systematically  eliminated  allowing  the  subsequent  exposure  of  weaker  signals 
and  their  emitters.  To  avoid  this,  the  communications  might  have  to  use 
bursts,  such  that  the  useful  integration  time  T  is  limited,  in  accordance  with 
increasing  the  number  of  emitters  in  the  network,  such  that 


-r  ^ 


(3.45) 
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This  also  implies  that  the  emitter  transmit  duty  cycle  would  need  to  be 
increased,  perhaps  using  dummy  traffic,  to  nearly  unity.  This  could  have 
operational  consequences  as  well  as  shorten  battery  recharge  intervals 
significantly. 

Since  the  output  of  the  cross  correlator  at  x  *  is  the  power  of 

m 

sm(t),  it  is  not  possible  to  avoid  the  cross  correlation  spike  without 

reducing  transmitted  power  once  antenna,  propagation  and  masking  issues  are 

fixed.  However,  it  may  be  possible  to  significantly  degrade  the  operation  of 

a  network  of  cross  correlators  by  using  repeated  short  codes  in  the  DS 

signal.  This  would  cause  a  number  of  minor  cross  correlation  spikes  with 

increasing  amplitude  until  the  spike  of  height  S  is  achieved  at  variable 

m 

t  *.  The  many  minor  spikes,  however,  could  make  it  very  difficult  to  pass 
reliable  emitter  Line  of  Bearing  (LOB)  information  between  intercept  receivers 
to  achieve  emitter  location.  Alternately,  all  emitters  could  share 
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a  certain  code  segment  or  segments  which  would  generate  many  false  spikes  and 
make  the  values  very  difficult  to  determine.  Sorting  on  only  the 

largest  spikes  would  allow  the  emitters  with  lower  power  to  remain  masked. 
Again,  however,  destruction  of  the  emitters  with  the  strongest  power  would 
have  to  be  avoided  by  perhaps  using  power  control  and  rapidly  moving  the 
emitters  with  the  largest  power. 

Ultimately,  however,  the  masking  techniques  which  most  capitalizes  on  the 

inherent  weakness  of  the  non-linear  recievers  is  to  locate  DS  signals  where 

many  interferers  also  operate.  Of  course,  the  environment  cannot  be  so  dense 

that  the  communicators  cannot  synchronize  and  exchange  data.  This  leads  to 

the  second  case  of  interest  in  which  the  total  distortion  power  at  the  delay 

x  *  is  increased  by  the  sum  of  D__,  Dv._  and  D__.  This  yields  the 
o  SP  NP  PP  1 

output  SNR  as  given  by  Equation  (3.41),  which  may  be  rewritten  in  the  form 
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where  is  defined  to  equal  the  spectral  taper  factor 


sin  x. 

2 

sin  ir  (f  -  f.  )T 

k 

c  k  c 
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k  J 

w  c  k  c  j 

where  f  is  the  center  frequency  of  the  narrowband  interferer  and  T  is 
the  chip  duration  of  the  DS  communications.  The  average  value  of  is 
approximately  1/2.22  which  will  be  used  to  simplify  Equation  (3.46). 


Note  from  Equation  (3.46)  that  the  interferers  have  a  dramatic  impact  on 
the  output  SNR.  For  instance,  neglecting  the  contributions  from  DgS»  Dgp 
and  D^p  yields  an  output  SNR  equal  to 
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(3.48) 


which  is  expected  to  be  low  for  reasonable  values  of  K  and  P 


k' 

With  Nj  much  less  than  the  typical  signal  power,  S^,  Equation  (3.46) 
becomes 
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(3.49) 


Clearly,  the  cross  correlation  receiver  must  try  to  eliminate  or  excise  the 
narrowband  interference  in  the  input  before  multiplication  or  after 
integration  or  both.  However,  in  a  dynamic  environment,  the  integration  time 
T  may  be  long  compared  to  the  average  time  interval  in  which  the  number  and 
total  power  of  the  interferers  changes  significantly.  This  makes  excision  at 
the  integrator  output  more  difficult. 


The  following  definitions  simplify  the  evaluation  of  Equation  (3.49) 


a  =  S  /S  (3.50a) 

m  mo 


and 


(3.50b) 


Using  Equation  (3.50),  Equation  (3.49)  becomes 
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(3.51) 


Generally,  if  the  output  SNR  for  a  DS  is  several  dB,  say  8  to  18,  then 
the  presence  of  the  autocorrelation  spike  can  be  detected  and  the  setting  of 
the  variable  delay  at  detection  provides  LOB  information.  As  mentioned 
earlier,  one  means  of  complicating  this  operation,  and  thus  degrading 
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performance,  is  to  create  compound  signals  with  multiple  autocorrelation  and 
cross  correlation  spikes. 


For  approximately  equal  powers  of  spread  and  unspread  (narrowband)  signals 


otm  equals  about  (3  which  equals  about  unity  for  all  values  of  m  and  k. 


In  this  case.  Equation  (3.51)  becomes 


1.47  (M2  -  M)  +  1.80  MK  +  KTW 


[1.47  +  1.80M] 

ss 


This  expression  confirms  that  narrowband  interference  needs  to  be  excised  from 
the  input  or  appropriately  filtered  or  ignored  at  the  output  of  a  cross 
correlation  receiver  in  order  for  it  to  achieve  reliable  performance. 


3.3.3  Detection  Using  a  Chip  Sate  Detector 


As  indicated  earlier,  the  potential  exists  for  the  cross  correlation 


receiver  to  develop  a  LOB  on  each  DS  emitter  for  which  the  differential 
propagation  between  antennas  is  more  than  about  (W  sfl ( SNR )  )  1 


seconds.  This  timing  difference  would  typically  correspond  to  a  few  feet  with 


W  equal  to  100  MHz  or  a  few  tens  of  feet  when  W  equals  10  MHz.  If  a 
s  s  th  ®  s 

feature  of  the  m  LOB  is  based,  then  multiple  receivers  and  feature 


extractors  could  isolate  individual  DS  emitters  unless  they  were  very  close 
together.  In  terms  of  targeting,  such  a  cluster  could  be  a  high  value 
target.  One  feature  which  might  possibly  be  used  as  a  tag  to  identify  a  LOB 
is  the  emitter  chip  rate.  However,  the  chip  rates  are  expected  to  be  very 
nearly  identical.  The  issue  then  is  how  well  can  a  chip  rate  detector 
distinguish  between  the  very  nearly  chip  rates  used  by  the  emitters  in  a  DS 
network.  This  issue  is  discussed  below  and  a  detailed  analysis  of  the 
performance  of  a  chip  rate  detector  is  provided  in  Appendix  A,  which  was 
prepared  by  Edgar  German  of  GT-Tech,  Inc. 


A  functional  block  diagram  of  a  chip  rate  detector  is  shown  in  Figure  3-9. 
The  basic  goal  of  the  receiver  is  to  multiply  the  receiver  input  by  a  delayed 
version  of  the  input  and  then  bandpass  filter.  Note  this  could  be  achieved 
using  two  channels  after  has  been  determined  by  placing  an  additional 

delay  in  one  channel  below  multiplication  and  bandpass  filtering. 

The  principle  of  operation  for  the  chip  rate  detector  is  as  follows. 

Given  a  chip  sequence  u  (t)  the  product  u  (t)u  (t-t)  is  random. 

o  o  o 

However,  the  expected  value  of  u  (t)u  (t-x)  is  a  deterministic 

oo 

square-wave  function  with  a  period  equal  to  the  chip  duration.  Since  the 

expected  value  of  u  (t)u  (t-t)  in  the  autocorrelation  function,  the 

0  0 

power  spectrum  of  the  product  has  harmonics  at  multiples  of 

the  reciprocal  of  the  interval.  The  fundamental  frequency  of  the  discrete 
power  spectrum  is  the  chipping  rate.  Hence,  extraction  of  this  spectral  line 
and  the  measurement  of  its  frequency  provides  a  measure  of  the  chip  rate 
transmitted  by  a  DS  emitter.  Stochastic  signals  which  have  periodic 
statistics,  such  as  the  autocorrelation  function,  are  called  cyclostationary 
signals . 

Three  factors  dominate  the  distinguishing  capability  of  the  chip  rate 
detector.  One  is  the  narrowband  interference  whose  dramatic  effect  was 
discussed  for  the  cross  correlator.  Narrowband  interference  has  similar 
effects  on  the  output  SNR  of  the  chip  rate  detector  as  it  did  for  the  cross 
correlator.  The  second  factor  is  the  filter  resolution  of  the  chip  rate 
detector.  The  processing  time  T  must  be  such  that  the  bandwidth  of  the 
bandpass  filter  1/T  is  less  than  or  equal  to  the  difference  is  the  chip  rates 
to  be  distinguished.  For  instance  a  rate  difference  of  one  chip  per  second 
requires  a  filter  with  one  Hertz  of  bandwidth  or  more.  The  third  factor  is 
the  extent  of  self-masking  interference  the  several  chip  sequences  have  on 
each  other.  These  last  two  factors  are  treated  in  detail  in  Appendix  A  and 
the  self-masking  factor  is  discussed  below. 


From  Appendix  A,  the  output  SNR  for  the  m  DS  signal  is  given  by 
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where  it  is  assumed  that  the  signal  power  of  the  DS  signals  is  well  above  the 

thermal  noise  level,  N  .  From  Equation  (3.53),  it  is  clear  that  with  only 

one  DS  signal  present  the  output  SNR  is  bounded  at  the  value  TW.  The  reason 
2 

why  the  term  Sq  appears  in  the  denominator  is  that  even  when  the  DS 

signal  is  very  strong  it  is  still  random  and  this  random  behavior  generates 

2 

self  noise  which  has  a  power  of  Sq/TW  in  the  output  bandwidth  of  1/T 
Hertz  centered  on  the  frequency  W  Hz. 


When  all  the  M+l  signal  powers  are  the  same,  the  output  SNR  is  the  same 
for  all  signals  at  the  value  of 
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Since  the  signal  of  interest  is  a  sinusoid  in  essentially  Gaussian  noise,  an 
envelop  detector  at  the  output  of  the  bandpass  filter  will  reliably  detect  the 
chip  rate  line  when  (SNR)q  is  greater  than  about  10  dB.  The  probabilities 
of  detection  and  false  alarm  are  conveniently  related  to  the  value  of  the 
output  SNR  in  Appendix  VII  of  the  report  included  here  as  Appendix  A.  For  the 
output  SNR  to  be  less  than  10,  Equation  (3.54)  shows  that  the  number  of 
masking  emitters,  M,  must  satisfy 


M  > 


(3.55) 


For  values  of  TW  greater  than  10  ,  the  number  of  masking  emitters  must 
exceed  817.  Hence,  it  would  appear  that  emitters  are  not  able  to  mask  each 
other  from  a  chip  rate  detector  unless  the  chip  rates  are  very  close  to  one 
another.  Certainly  the  rates  should  differ  by  no  more  than  about  1/10  of  a 
chip  per  second.  This  means  that  for  moving  emitters  or  receiver  the  time 
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dilation,  the  source  of  doppler  shift,  will  often  provide  more  differences  in 
the  apparent  chip  rate  than  the  transmitter  timing  sources.  Such  timing 
sources  are  extremely  stable.  If  the  chip  rates  differ  by  more  than  a 
fraction  of  a  chip  per  second,  then  filtering  may  be  able  to  isolate  the 
individual  chip  rate  lines  in  bandwidths  where  the  SNR  is  quite  large  for  even 
relatively  large  networks.  As  in  the  case  of  the  cross  correlator,  it  appears 
that  narrowband  interference  is  the  key  to  whether  the  chip  rate  detector  will 
achieve  the  necessary  values  of  output  SNR. 

One  possibility  to  lower  the  risk  of  being  isolated  by  a  chip  rate 
detector  is  to  jitter  the  chip  rate.  However  the  amount  of  jittered  required 
is  significant  and  therefore  the  structure  of  the  despreading  correlation  in 
the  friendly  receiver  is  complex,  unless  a  convolver  implementation  is  used. 

A  second  possibility  is  to  limit  the  integration  time  T  to  a  maximum  usable 
value  of  about  1/2R  by  hopping  at  a  rate  of  R  hops  per  second.  This  has 
significant  potential  as  just  a  few  tens  of  hops  per  second  will  dramatically 
reduce  the  usable  integration  times  or  increase  the  usable  filter  bandwidths. 
Hence,  a  hybrid  or  jittered  DS  system  or  both  should  be  considered  when  the 
goal  is  to  avoid  a  persistent,  high  technology,  non-linear  intercept  receiver 
threat. 

Alternately,  the  chip  shape  may  be  selected  to  radically  reduce  or 
eliminate  the  chip  rate  line.  Theoretically  this  happens  for  sin  x/x 
time-shaped  pulser,  which  have  a  uniform  bandlimited  Fourier  spectrum^2^. 

While  those  pulses  are  not  realizable,  the  theoretical  null  for  uniform 
spectrum  can  be  approached  with  realizable  pulse  shapes.  In  this  way,  the 
chip  rate  line  can  be  very  much  reduced.  As  a  practical  matter,  chip  shapes 
which  are  bandlimited  have  excessive  temporal  content  at  times  well  outside 
the  "chip"  interval.  However,  compromises  and  partial  response  techniques  may 
exist  allowing  significant  reduction  in  the  performance  of  a  chip  rate 
detector. 


[2]  Reed,  D.E.,  and  Wickert,  M.A.,  "Minimization  of  Detection  of  Symbol-Rate 
Spectral  lines  by  Delay  and  Multiply  Receivers,"  IEEE  Trans,  on  Comm.,  Vol. 
36,  No.l,  January,  1988,  pp.  118-120. 
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Of  course,  placing  the  DS  signals  in  a  dense  signal  environment  will 
dramatically  degrade  the  performance  of  the  chip  rate  detector  and  the  cross 
correlation  receiver.  In  this  case,  however,  some  of  the  "TW"  product  of  DS 
signal  must  be  used  to  overcome  the  power  of  the  in-band  interference.  Hence 
the  signal  power  spectral  density  can  not  be  as  low  compared  to  the  thermal 
noise. 


ssssss* 

3 

% [ 
1 

CHAPTER  4 


CONCLUSIONS 


The  technical  analyses  of  the  report  began  with  consideration  of  angle 
sorting  or  isolation  using  directional  antennas.  The  conclusion  of  the 
antenna  analyses  is  that  at  400  MHz  it  is  very  unlikely  that  intercept 
antennas  will  be  large  enough  to  resolve  individual  emitters.  In  fact,  from 
an  intercept  range  of  about  15  miles,  emitters  within  2  miles  of  one  another 
are  likely  to  be  unresolved.  For  a  given  aperture  size,  the  beamwidth  of  an 
antenna  at  1800  MHz  is  about  0.22  times  narrower  than  at  400  MHz.  Hence,  for 
a  15  mile  intercept  range,  the  resolution  is  about  0.44  miles.  These  numbers 
assume  a  25  foot  azimuth  aperture  which  is  large  for  an  airborne  antenna 
anticipated  when  the  intercepted  range  is  15  miles  or  more.  Airborne  arrays 
using  multiple  aircraft  are  a  potential  means  of  achieving  the  required 
effective  aperture  and  this  intercept  threat  should  be  considered  if  there  is 
any  evidence  of  its  existence. 

The  second  part  of  the  antenna  analysis  reveals  that  with  imminently 

usable  antenna  sizes  the  intercepted  power  for  DS  signals  may  be  well  above 

2 

the  thermal  noise  floor  for  free-space,  1/R  propagation.  Hence  while 

antennas  large  enough  to  resolve  emitters  in  angle  are  unrealistically  large, 

2 

antennas  of  realistic  dimensions  provide  enough  gain  over  1/R  paths  to 

4 

achieve  large  values  of  Signal-to-Noise  Ratio  (SNR).  However,  over  1/R 
propagation  typical  of  UHF  ground-to-ground  intercept  paths,  the  antennas 
again  becomes  very  large.  Being  on  the  ground,  however,  such  antennas  could 
be  built.  But  it  appears,  that  as  intercept  range  increases  to  protect  the 
receiver  antenna  from  hostilities  near  the  tactical  emitter,  the  required 
antenna  size  remains  too  large  to  survive  in  wartime.  Collection  sites  may 

4 

exist  m  peacetime  which  pose  threats  to  DS  signals  over  1/R  paths. 

The  typical  intercept  receiver  for  tactical  DS  signals  is  expected  to  have 
many  DS  signals  present,  each  one  being  large  compared  to  the  thermal  noise 
but  comparable  in  power  to  each  other.  The  intercept  receiver,  likely  to  be 
airborne,  is  also  expected  to  find  a  dense  environment  of  narrowband  signals 
which  result  in  a  significant  level  of  interference  power.  It  is  in  this 


intercept  context  that  the  majority  of  the  effort  documented  in  this  report 
was  focused. 

In  cnis  intercept  context  it  is  clear  that  the  power  of  any  one  signal  is 
typically  well  below  the  sum  of  the  powers  of  the  several  other  signals 
present  in  the  receiver  as  a  result  of  their  emitters  lying  within  the 
mainbeam  of  the  intercept  antenna.  Hence,  Fourier  analysis  or  linear 
filtering  by  any  method  are  not  expected  to  perform  well.  The  emitters  are 
said  to  mask  each  other  as  the  expected  intercept  SNR  is  on  the  order  of  1/M 
where  M  is  the  total  number  of  DS  emitters  in  the  antenna  beamwidth. 

Therefore,  non-linear  receivers  are  required  which  have  the  ability  to  perform 
when  the  power  of  the  Signal  of  Interest  (SOI)  is  well  below  the  sum  of  the 
powers  of  the  Signals  Not  of  Interest  (SNOI).  Conclusions  regarding 
performance  of  such  receivers  are  given  below. 

Non-linear  receivers  are  able  to  detect  periodic  and  certain  other  aspects 
of  a  signal  even  when  the  signal  power  is  less  than  the  noise  power.  This 
includes  the  energy  of  the  signal  at  zero  or  twice  the  frequency,  the  Angle  of 
Arrival  (AOA)  of  the  signal,  the  carrier  frequency  by  measurements  at  twice 
the  carrier  frequency  for  BPSK  spreading,  the  chip  rate  or  its  harmonies,  and 
the  hop  rate  or  its  harmonics.  Baseband  symbol  rates  can  also  be  useful  in 
distinguishing  between  the  SOI  and  SNOI.  Non-linear  receivers  include  the 
total  power  radiometer,  the  frequency  doubler,  the  auto  and  cross  correlation 
receivers,  and  a  variety  of  receiver  subsystems  called  feature  detectors. 

Since  non-linear  receivers  extract  signal  power  their  performance  is 
expected  to  be  good  when  there  is  only  one  signal  plus  thermal  noise. 
Furthermore,  the  cross  correlation  receiver  is  able  to  adjust  the  relative 
signal  delay  such  that  a  spike  corresponding  to  the  power  of  each  signal  may 
be  extracted  and  perhaps  isolated.  Emitters  whose  differential  delay  between 
the  antennas  of  a  dual  channel  cross  correlator  is  a  nanosecond  (foot)  or  more 
may  be  resolved  and  a  Line  of  Bearing  (LOB)  determined.  Hence  the  number  of 
DS  emitters  in  the  mainbeam  may  be  known  and  hence  emitter  masking  is  not 
complete.  Without  a  dense  signal  environment  or  attention  to  signal  design, 
emitters  within  networks  are  not  expected  to  be  closely  spaced  enough  that  the 


self  masking  of  DS  emitters  is  sufficient  to  overcome  the  ability  of  a  cross 
correlator  to  isolate  them.  Repeated  short  codes  for  a  given  emitter  and 
common  code  segments  among  emitters  produce  a  number  of  off-peak 
autocorrelation  spikes  and  cross  correlation  spikes.  While  the  largest 
autocorrelation  spikes  are  still  as  large  with  sub-coding  as  without,  the  many 
off-peak  spikes  greatly  complicate  the  detection  process.  More  importantly, 
the  correlation  functions  with  many  spikes  complicate  the  fundamental  problem 
of  tagging  LOB's  so  they  can  be  coordinated  between  receiver  sites  to  provide 
emitter  location.  Hence,  it  would  not  be  unexpected  for  a  dedicated,  highly 
computational  intensive  interceptor  operating  in  a  noise-only  environment  to 
derive  a  LOB  for  each  DS  emitter  in  a  network  of  DS  emitters.  However,  signal 
design  may  vastly  complicate  the  problem  of  usefully  coordinating  this 
information  between  receiver  sites  by  assigning  identifying  feature  tags  to 
each  LOB.  In  fact,  further  signal  design,  can  dramatically  reduce  the  threat 
that  useful  features  can  be  extracted  at  all. 

While  LOB  detection  can  be  made  complicated  through  design  of  signals  with 
command  and  composite  codes,  dramatic  loss  of  receiver  capability  is  not 
expected  until  further  degradation  is  caused  by  a  dense  environment  of 
unwanted,  narrowband  SNOI.  The  combination  of  signal  design  to  avoid  single 
spike  autocorrelation  functions,  design  to  produce  spikes  in  the  cross 
correlation  functions,  design  to  avoid  feature  extraction,  and  the  presence  of 
a  dense  environment  of  SNOI  should  provide  protection  of  the  emitters  in  a 
very  wide  variety  of  tactical  and  surveillance  situations. 

A  key  element  in  the  design  of  Low  Probability  of  Intercept  (LPI)  signals 
is  to  avoid  putting  features  into  signals  which  can  be  extracted  using 
nonlinear  methods.  For  instance,  the  performance  of  chip  rate  detectors  can 
be  limited  by  jittering  the  chip  rate  and  limited  much  more  by  Frequency 


ft 


Hopping  (FH)  since  this  limits  the  useful  bandwidth  of  the  bandpass  filter 
searching  for  the  chip  rate  "line.”  The  chip  rate  "line"  is  broadened  by  the 
hopping  to  a  null-to-null  bandwidth  equal  to  twice  the  hopping  rate. 
Furthermore,  the  Fourier  spectrum  of  the  chip  can  be  tailored  to  lower  the 
harmonic  content  at  the  chip  rate  or  any  other  frequency.  The  more 
bandlimited  the  chip  shape  is  the  less  harmonic  content  at  the  chip  rate. 

Such  signals,  however,  do  have  "time-tails"  which  complicate  friendly 
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reception  and  increase  the  drop  in  network  performance  and  capacity  due  to 
multipath  and  intersymbol  interference.  Again,  a  dense  environment  of 
narrowband  SNOI  has  a  dramatic  and  significant  impact  on  the  complexity, 
cornputer/processing  content,  and  performance  of  the  non-linear  chip  rate 
detector.  Without  specific  design  against  chip  rate  detectors  and  in  the 
absence  of  narrowband  interference,  it  appears  that  a  bank  of  bandpass  filters 
with  individual  bandwidths  of  W  Hertz  could  differentiate  and  isolate  emitters 
whose  chipping  rates  differ  by  W  chips/second.  There  are  no  "sum  and 
difference"  chip  rate  lines  as  might  be  initially  expected. 

Another  key  emitter  feature  is  the  exact  center  frequency  of  an  emitter. 

If  this  frequency  can  be  repeatably  differentiated  among  emitters  it  can  be 
used  as  an  LOB  tag  like  the  chip  rate  could.  However,  signal  design  can  be 
used  to  complicate  the  task  of  useful  carrier  recovery  by  a  noncooperative 
receiver.  Such  methods  include  nonuniform  envelope  and  uniform  envelope 
spreading  sequences  derived  by  digitizing  bandlimited  noise,  various  forms  of 
QPSK  spreading  FH  to  avoid  narrowband  filtering  (long  integration  time).  Time 
Hopping  (TH)  to  avoid  long  integration  times,  and  jitter  of  the  center  carrier. 

As  an  overall  result,  it  appears  that  against  conventional  Fourier 
receivers  that  DS  emitters  within  a  DS  network  do  in  fact  mask  each  other. 
Against  a  cross  correlation  or  cross  power  spectrum  analyzer  the  emitters  with 
the  beamwidth  may  be  counted.  Common  and  compound  short  codes  complicate  the 
required  detection  and  a  dense  electromagnetic  environment  plus  appropriate 
signal  design  should  make  LOB’s  very  difficult  to  extract  except  for  emitters 
with  very  large  received  SNR.  Tagging  these  LOB's  with  useful  and  reliable 
features  can  be  made  very  unlikely  when  signal  design  procedures  are  used  and 
the  operational  environment  is  dense  with  dynamically  varying  SNOI. 
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ANALYSIS  AND  SIMULATION  OF  A  CHIP  RATE  DETECTOR 
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WITH  A 

MULTIPLICITY  OF  PSEUDONOISE  INPUTS 
1.0  INTRODUCTION 


A  standard  method  of  detecting  the  presence  of  a  pseudorandom  (PN), 
spread  spectrum,  waveform  is  through  the  use  of  a  chip  rate  detector. 
However,  multiple  chip  stream  inputs  at  slightly  differing  chip  rates  may 
conceal  each  other,  or  at  best,  generate  "masking"  noise  for  each  other. 
The  analyses,  and  the  simulation  performed,  address  and  resolve  the 
following  issues: 

1.  Effects  of  multiple  chip  streams  upon  each  other  in  terms  of  the 
signal-to-noise  ratios  of  the  various  streams  in  the  detector 
output . 

2.  Resolution  between  chip  streams  at  different  rates,  but  with  the 
same  encoding. 

3.  Masking  effects  of  multiple  chip  stream  interference  upon  the 
detectability  of  an  individual  chip  stream. 

4.  Effect  of  receiver  noise  upon  the  detectability  of  multiple  chip 
streams. 

Simulation  has  been  used  to  support  the  analyses  by  verifying  the 
analytical  results.  Analyses  to  resolve  these  issues  has  involved: 

1.  Identification  of  the  chip  rate  spectral  lines  of  the  detector 
output  for  a  single  chip  stream  as  a  reference  case. 

2.  Identification  of  the  continuous  spectrum  of  the  detector  output 
for  a  single  chip  stream  as  a  reference  case. 
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3.  Identification  of  the  chip  rate  spectral  lines  of  the  detector 
output  for  multiple  chip  streams,  including  any  "mixing"  harmonics 
(if  present). 

4.  Identification  of  the  continuous  spectra  of  the  detector  output 
for  multiple  chip  streams,  including  any  "mixing"  harmonics. 

5.  Calculation  of  the  ratios  of  the  continuous  and  discrete  spectral 
components . 

6.  Calculation  of  the  output  signal-to-interference  ratio  of  the 
detector  output  filter  as  a  function  of  the  bandwidth  (integration 
time)  of  that  filter.  This  calculation  considers  both  multiple 
chip  streams  and  receiver  noise  as  components  of  the  interference. 
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Modeling  of  the  chip  rate  detector  has  been  performed  at  baseband  to 
simplify  the  analysis.  In  block  diagram  form,  the  model  appears  in  Figure 
1-1 . 
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Figure  1-1  Chip  Rate  Detector 

Derivation  of  the  chip  rate  detector  output  characteristics  will 
involve  computation  of  the  multiplier  output  ensemble  mean  frequency  content 
and  power  spectrum.  To  this  end  the  mean  and  autocorrelation  function  of 
the  output  will  be  computed  and  Fourier  transformed  as  a  simpler  means  of 
calculation. 


2.0  SUMMARY  OF  RESULTS 


While  the  analyses  to  determine  the  chip  rate  detector  output 
signal-to-noise  ratio(s)  are  quite  lengthy,  the  basic  results  are  simple  to 
write  and  to  interpret.  Sections  1  to  5  describe  the  results  and 
interpretation  of  the  investigations  of  those  results  and  form  a 
self-standing  document.  Section  6  and  the  Appendices  detail  the  derivation 
of  the  basic  results. 

Output  spectra  (chip  rate  lines  and  noise)  of  the  chip  rate  detector 
are  of  three  basic  types: 

I .  Individual  Spectral  Components 

Each  chip  stream,  or  noise,  has  a  pair  of  components  in  the  detector 
output  that  are  the  components  that  would  be  present  for  that  individual 
input  alone.  Components  for  the  individual  inputs  are: 

Chip  Streams 

a.  A  periodic  square  wave  at  the  component  chip  rate,  which  produces 
the  desired  chip  rate  spectral  lines. 

b.  A  continuous  spectrum  that  is  uncorrelated  with  the  periodic 
square  wave,  and  has  a  width  of  twice  the  chip  rate  at  the  first 
nulls.  This  spectral  component  causes  "self"  noise  in  the 
detector  output. 

There  are  no  mixing  chip  rate  lines:  lines  that  could  possibly  occur 

due  to,  say,  the  sum  and  difference  of  the  different  chip  rates. 

Noise 


a.  An  uninteresting  DC  component  that  is  the  value  of  the  noise 

autocorrelation  function  at  the  half  chip  delay. 

b.  A  continuous  spectrum  that  is  basically  the  convolution  of  the 

noise  spectrum  with  itself.  It  has  a  width  of  approximately  twice 
the  chip  rate.  Two  input  spectra  were  used  in  the  analyses:  a 
rectangular  spectrum  and  a  sin(x)/x  spectrum.  Use  of  simulation 
verified  the  results  for  both  noise  spectra.  In  practice,  the 

chip  rate  rectangular  spectrum  is  the  recommended  filter  as  it 
will  not  reduce  the  desired  chip  rate  lines  as  much.  For 

simplicity  in  the  analytical  calculations,  the  sin(x)/x  spectra 
was  chosen  for  the  complete  output  spectrum  calculations.  But,  at 
the  important  chip  rate  line,  the  rectangular  filter  noise 
contribution  was  also  calculated  exactly. 

For  N  chip  streams  at  differing  chip  rates,  there  are  N  discrete  line 
sets  at  the  various  chip  rates  and  N  sets  of  self  noise,  plus  the  noise  DC 
and  continuous  spectrum.  Assuming  the  chip  rates  of  all  the  chip  streams 
are  similar,  a  narrowband  output  filter  that  is  tuned  to  a  particular  chip 
rate  will  have  only  that  chip  rate  line  as  its  output  (given  sufficient 
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resolution).  All  the  continuous  spectral  components  (within  the  filter 
bandwidth)  will  also  appear  as  "noise"  in  the  filter  output. 


Resolution  between  two  different  chip  rate  lines  requires  a  filter 
whose  bandwidth  is  less  than  the  separation  between  the  chip  rates. 
Additionally,  the  filter  must  have  enough  of  a  roll-off  to  reject 
neighboring  chip  lines  over  a  meaningful  dynamic  range. 


II .  Intermodulation  Components 


Each  chip  stream  interacts  with  every  other  chip  stream  and  the  noise. 
These  interactions  occur  pairwise.  They  are  approximately  twice  the  chip 
rate  in  bandwidth  and  appear  as  "noise"  in  the  narrowband  filter  output. 


The  number  of  pairwise  components  are  the  the  number  of  pairs  that  can 
be  formed  from  N  +  1  (N  chip  streams  plus  noise)  things: 


(chip/chip)  +  N  (chip/noise) 


NfN+11 


which  increases  as  the  square  of  the  number  of  chip  streams, 


III.  DC  Intermodulation  Components 


DC  components  representing  the  energy  of  the  interchip  modulation  are 
present  in  the  output  of  the  multiplier.  They  are  uninteresting  as  they  do 
not  appear  in  the  chip  rate  detector  output  filter. 


and  Noi; 


iignal-to-Nois 


By  collecting  the  aforementioned  terms  (Section  6  and  the  Appendices) 
the  chip  rate  detector  output  signal-to-noise  ratio  can  be  calculated. 
Basic  assumptions  for  the  relationship  are: 


1.  The  narrowband  filter  has  a  noise-bandwidth  B  =  1/T  that  is 
sufficiently  narrow  to  reject  other  (different)  chip  rate  lines. 


2.  Chip  rates  of  the  different  streams  are  close  to  one  another  (say, 
within  107.)  so  that  output  noise  spectra  are  relatively  the  same 
for  pairwise  intermodulation  products. 


The  signal-to-noise  ratio  for  a  particular  chip  stream  output  with  a  pair  of 
chip  stream  inputs  and  noise  is: 
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with 

S^2^  =  input  power  of  chip  stream  1(2), 
and 

N-  =  input  noise  power. 

and  is  the  basic  result  of  this  report! 

The  constants  a  and  P  are  highly  dependent  upon  the  modeling  of  the 
noise  spectrum  used.  A  low-pass  rectangular  input  filter,  prior  to  the  chip 
rate  detector,  with  a  bandwidth  equal  to  the  chip  rate  requires  a  =  1,  and 
agrees  with  textbook  analyses  for  a  single  chip  stream  and  noise.  More 
conveniently,  for  analytical  purposes  and  as  another  case  for  verifying  the 
basic  methodology  by  simulation,  is  filtering  of  only  the  noise  with  a 

9 

filter  having  the  chip  shape  as  an  impulse  response  requires  a  =  6/ir'  . 
This  choice  for  a  results  in 

7}  3 

a  T  =  7’ 

and 

3  =  1. 

In  tabular  form: 


Both  noise  spectra  were  included  in  the  simulation  and  verified  within  1  dB. 


Results  can  also  be  written  in  terms  of  signal-to-noise  ratios  of  the 
input  signals,  but  interpretation  in  that  form  will  be  deferred  for  the 
moment.  The  other  chip  stream  has  the  same  expression  for  the  output 


signal-to-noise  ratio;  but.  with  the  numorat.or  identifying  that  rhip 
amplitude  (denominators  are  identical). 

For  interpretation  of  the  result,  two  cases  are  of  interest;  large  and 
small  input  signal-to-noise  ratios. 

2.1.1  Small  Input  Sienal-to-Noise  Ratio 

Given  the  input  signal-to-noise  ratio  for  the  reference  chip  stream 


then: 
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and  the  result  depends  upon  the  ratio  of  the  two  chip  stream  powers. 
Case  I:  The  Chip  Stream  S^  Dominates: 


S1  >>  S2 
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which  is  the  well  known  result  for  a  chip  rate  detector  with  only  a  single 
chip  stream  and  noise  in  the  input  (a  =  1). 

Case  II:  The  Chip  Stream  S9  Dominates: 
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where: 


S' 

Mi,  2 


A  S2 


and,  the  results  are  written  in  terms  of  the  signal-to-noise  ratios  of  the 
input  chip  streams. 


Two  subcases  present  themselves  according  to  whether  the 
signal-to-noise  of  the  second  chip  stream  is  small  or  large.  If  the  second 
chip  stream  has  a  small  signal-to-noise  then  the  output  signal-to-noise 
ratio  reduces  to  simply  the  single  chip  stream  case  as  expected. 

For  a  large  signal-to-noise  of  the  second  chip  stream,  the  output 
signal-to-noise  ratio  becomes: 


2 

which  is  worse  than  simply  noise  by  the  factor  of  a  x  /4. 
2.1. 2  Large  Input  Signai-to-Noise  Ratio 


i'- 
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For  a  large  input  signal-to-noise  ratio 


»  1, 

rewrite  the  expression  above  as: 
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1  +  Sg/S^j  +  3S2/S1  +  S2/S1"/J[N/S]i "i’ 


which  depends  upon  the  ratio  of  the  power  in  the  two  chip  streams.  Again, 
two  cases  present  themselves. 

Case  I:  The  Chip  Stream  Dominates: 

In  this  instance  the  apparent  output  is  that  of  the  reference  chip 
stream  by  itself.  Only  the  self-noise  of  the  reference  chip  stream  is 
present.  And, 

S1  >;>  S2 
and 


o,l 


=  TV, 


The  output  is  simply  the  "processing  gain"  of  the  detector. 
Case  II:  The  Chip  Stream  S2  Dominates: 
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With 


S2  >:>  S1 


S' 

y  o,l 


The  output  signa’l-to-noise  ratio  depends  only  upon  the  ratio  of  the  chip 
stream  powers,  as  might  be  expected. 

2.2  Generalization  to  N  Chip  Stream  Inputs 

Calculation  of  the  signal  and  noise  power  for  N  chip  streams  in  the 
chip  rate  detector  input  is  a  relatively  simple  generalization  of  the  two 
chip  stream  case.  From  Appendix  VI,  the  total  "noise"  out  of  the  chip  rate 
detector  is: 


■I  ]  *  * r'2/4  NM 1 1  Vk +  3^X  sjni. 


j  k 


where  the  sums  are  over  the  N  chip  streams. 

By  making  reasonable  assumptions  about  the  distribution  of  the  powers 
(signal-to-noise  ratios)  of  the  chip  stream  inputs  an  average  output  noise 
can  be  calculated.  This  average  is: 

E[TN]  =  N2  N  E  [SNR2]  +  a  tt2/4  +  |  N(N-l)  E  [SNR] 2  +  3/3N  E  [SNR] 


The  output  chip  line  power  (normalized  to  the  input  noise  power)  for  the  j 
chip  stream  is 

S  =  S2/N2  TV 

with  the  average  output  equal  to: 

E[S]  =  E[SNR2]  . 
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With  these  quantities,  the  average  output  signal-to-noise  ratio,  over 
all  chip  streams  is: 


SNR 


av 


TW 

1  +  a  x2/ (4N)  +  |  ( N-l )  E[SNR]2/E[SNR2]  +  3^E [SNR] /E [SNR2] 


Further  evaluation  of  the  average  signal-to-noise  ratio  requires 
specification  of  the  distribution  function  of  the  chip  stream 
signal-to-noise  ratios. 

Two. cases  will  be  considered.  One  will  be  a  uniform  distribution  over 
dB  of  the  signal-to-noise  ratios  of  the  chip  stream  transmitters.  The  other 
will  be  a  uniform  distribution  over  transmitter  ranges,  with  all  transmitted 
powers  equal. 

Case  I : 

Cniform  Distribution  of  Signal-to-N'oise  Ratios  fin  dB) . 

Assume  the  distribution  is  given  between  zero  and  dBmax,  then 


pdB(dB)  =  an 


max 


fl,  0  <  dB  <  dB 
0,  otherwise. 


max’ 


and  obtaining  the  distribution  of  S  through  the  relation 


dB  =  10  log1Q(S): 


PS^  =  in  5  5’  1  -  S  -  S(nax- 

max 

With  these  relations,  the  average  chip  rate  detector  output  signal-to-noise 

ratio  is  (for  S  „  >>  1) : 

'  max  ' 


SNR 


_  1 
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1  *  a  t2/(4N)  *  3 (N— I ) / In  SBax  *  6/J/SBax 


For  a  large  number  of  emitters  (N  >>  1)  and  high  input  maximum 
signal-to-noise  ratios  (Smo  >>  100  (20  d  B)),  the  the  formula  may  be 

approximated  by: 
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1  +  3 (N— 1 ) /In  Sn 


A  conclusion  can  be  reached.  The  average  signal-to-noise  ratio 
decreases  at  least  as  fast  as  N  for  emitters  spread  over  extremely  wide 
dynamic  ranges.  For  a  60  dB  dynamic  range: 


lav  =  ?T 


1  +  .22  (N-l) 


Case  II: 

Uniform  Distribution  of  Rant 


Assume  the  distribution  is  given  between  R  •  and  R _ ,  then 

°  mm  max 

, HFn  1  f1,  Rmin  -  R  -  Rmax’ 

p.n(dB)  =  n - 37-TT - 

'"ax  min  |0,  otherwise, 

and  obtain  the  distribution  of  S  through  the  relation 

s  =  4- 

As  in  case  I,  for  large  max  to  min  range  range  ratios: 


;av  =  !¥ 


1  *  a  *1 2/(4N)  ♦  f  (N-l)  R„in/Rmax*  W/SNR^ 


mm  range 


A»ain,  assuming  the  minimum  range  signal-to-noise  ratio  is  large,  and  N 
also: 


1  +  7  (N-l)  R  •  /R 

2  ’  min'  max 


^  «i  ■  w*  w  w~a 


A  similar  conclusion  to  the  uniform  distribution  in  dB  case  can  be 
reached.  The  average  signal-to-noise  ratio  decreases  as  the  number  of 
emitters  (N).  A  choice  of  maximum-to-minimum  range  of  1000  (60  dB 

dispersion)  gives,  in  this  instance: 


SNR 


TV 


av  "  S 


1  +  ,0045(N-1) 


Conclusion 

A  large  number  of  emitters  tends  to  obscure  each  emitter's  chip  rate 
line.  The  obscuration  increases  as  10  log(N)  in  dB  when  the  emitter 
signal-to-noise  ratios  are  distributed  over  a  wide  range  of  values.! 
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3.0  BASICS  OF  THE  ANALYTICAL  METHOD 


Analytical  techniques  used  within  this  report  derive  two  basic 
quantities;  even  though  the  derivations  are  extensive.  These  two  quantities 
are  the  amplitude  of  the  chip  rate  lines  and  the  noise  variance  (both 
receiver  noise  and  interchip  mixing  noise)  out  of  a  bandpass  filter  placed 
about  the  chip  rate  lines.  Both  quantities  are  obtained  by  calculating  the 
time  averaged  autocorrelation  function  and  its  Fourier  transform,  the 
spectrum,  of  the  chip  rate  detector  output.  The  autocorrelation  is  the 
primary  quantity  calculated,  as  it  reveals  the  cyclostationarity  of 
pseudonoise  chip  streams  (and,  therefore,  the  chip  rate  lines)  much  more 
readily  than  a  direct  (but  equivalent)  spectrum  calculation. 

These  quantities  are  not  random  variables;  they  represent  the  mean 
(amplitude)  of  the  chip  rate  sine  wave  and  the  noise  standard  deviation. 
Assuming  a  Gaussian  probability  density  for  the  total  noise,  the  probability 
of  detection  of  the  chip  stream  can  be  calculated  (See  Section  4.). 

The  spectrum  of  the  chip  rate  detector  output  (computed  via  the  Fourier 
transform  of  the  time  averaged  autocorrelation  function)  allows  the  signal 
amplitude  and  noise  variance  to  be  calculated  by  virtue  of  integrating  the 
spectrum  over  the  frequency  range  of  interest. 

Given  an  output  spectrum  S  (u;)  (a  noise  power  density  which  will  be 

Z 

calculated)  of  the  delay  and  multiply  circuit,  the  output  power  of  the 
narrowband  filter  H(w)  (of  bandwidth  B)  is: 


N 

o 


“00 


Sz(w)  d u 


If  the  filter  is  chosen  to  be  rectangular  and  have  a  narrow  bandwidth  of  B, 
about  some  frequency  uiQ,  then  the  output  noise  power  is: 


+  the  negative  frequency  contribution. 

As  both  chip  rate  sine  waves  and  noise  have  the  same  contributions  for 
negative  frequencies  (spectra  are  even  functions),  the  negative  frequency 
contributions  will  be  ignored  hereafter  as  only  ratios  are  of  interest. 

Assuming  that  the  spectrum  is  reasonably  smooth  within  the  confines  of 
the  filter,  then  the  noise  power  output  is: 


N 

o 
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Now,  B  is  in  radians  per  second.  Rewriting  the  expression  in  Hertz 

(B  =  B/2t)  : 

vo  = 8  SZK>  ■ 

=  T  W  ’ 

as  T  =  1/B  is  the  integration  time  of  the  filter. 

If  the  spectrum  includes  a  delta  function  at  the  frequency  uQ,  then 

the  filter  output  will  also  include  a  constant  signal  term  independent  of 
the  filter  bandwidth  (as  long  as  the  signal  term  falls  anywhere  within  the 
filter  bandwidth).  Let  the  sine  wave  spectrum  be  represented  by: 

s\(u0)  =  2*  A  6(u)  -  u>Q), 


then,  the  output  is 
SA  =  A, 


if  w'Q  -  B/2  <  u  <  uQ  +  B/2  zero  otherwise. 

These  relationships  are  used  to  interpret  the  results  for  the  spectra 
computed  subsequently. 

Actual  Filtering 

While  the  perfectly  rectangular  filter  gives  a  quick  means  of 
interpreting  the  results,  actual  filters  differ  in  detail.  The  simulation, 
used  in  developing  this  report,  uses  a  Fourier  transform  as  the  output 
filter  bank.  A  straight-forward  Fourier  transform  over  N  chips  is  identical 
to  using  an  integrate  and  dump  circuit  with  an  integration  time  of  NTc>  Or, 

as  M  consecutive  transforms  are  taken  with  their  magnitude  squared  output 
average  over  the  M  transforms,  the  integrate  and  dump  can  be  looked  upon  as 
a  sliding  window. 

The  noise  bandwidth  of  the  sliding  window  integrator  is  exactly  equal 
to  B  =  1/T,  the  ideal  filter  bandwidth.  However,  the  discrete  Fourier 
transform  acts  (with  the  equivalent  sliding  window)  as  though  it  has  a 
sin(x)/x  transfer  function  with  considerable  leakage  into  adjacent  channels 
(Fourier  coefficients).  It  is  also  difficult  to  read  between  the  "cracks" 
when  an  input  sine  wave  falls  between  the  discrete  transform  frequencies. 

By  using  a  bell-shaped  weighting  over  each  span  of  chip  detector  output 
of  T  seconds,  the  effect  is  a  bell-shaped  impulse  response  of  T  seconds 
duration.  The  well-known  cosine-squared  weighting  was  used  for  the 
weighting.  It  effectively  increases  the  noise  bandwidth  over  the  constant 
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amplitude  window  bv  a  factor  of  3/2.  The  signal-to-noise  ratio  of  each 
transform  channel  (filter)  is  decreased  by  3/2.  Leakage  of  the  adjacent 
channels  is  still  present;  but  the  outputs  are  more  readable  due  to  the 
wider  bandwidth  and  lack  of  perturbation  from  distant  channels. 
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4.0  PRnOAB  I LITY  (IF  DETECTION  OF  THE  CHIP  RATE  DETECTOR  PIT  PIT 


Efforts  described  in  this  report  are  devoted  to  computation,  and 
verification,  of  the  basic  signal-to-noise  output  of  the  chip  rate  detector 
under  various  conditions.  This  quantity  is  not  a  random  variable.  It  is 
the  ratio  of  the  mean  output  (squared)  to  the  noise  variance.  These 
quantities  specify  two  moments  of  the  probability  density  of  the  chip  rate 
detector  output.  If  the  output  is  assumed  to  be  the  sum  of  a  sine  wave  plus 
a  Gaussian  random  variable  (the  noise)  then  the  signal-to-noise  ratio 
calculated  in  this  report  specifies  the  random  process  in  its  entirety. 

As  the  noise  (chip  mixing  plus  receiver  noise)  is  narrowband  filtered; 
the  Central  Limit  Theorem  allows  us  to  infer  that  the  chip  rate  detector 
output  is,  in  fact,  composed  of  Gaussian  noise.  The  correlated  delta 
function  is  the  sine  wave  output.  An  envelope  detector  (placed  at  the  chip 
rate  detector)  output  will  be  the  optimum  estimator  of  the  sine  wave 
presence. 

Each  integration  time  of  the  bandpass  filter  (T)  produces  a 
statistically  independent  sample  of  the  sine  wave  plus  noise  process.  By 
placing  a  threshold  at  the  envelope  detector  output  the  presence  (or 
absence)  of  a  chip  line  can  be  estimated.  Two  quantities  express  the 
effectiveness  of  the  detector:  (1.)  the  probability  of  detection  (P^) ,  and 

(2.)  the  probability  of  a  false  alarm  (j3) .  These  quantities  are  derived 
from  Marcum's  Q  function  using  the  signal-to-noise  ratio  out  of  the  bandpass 
filter  as  the  argument.  A  nomogram  is  given  in  Appendix  VII  for  computing 
the  probabilities,  given  the  signal-to-noise  ratio;  or  any  quantity  given 
rhe  other  two.  As  an  example  from  the  nomogram: 


J  =  10  2 , 

then 

SNR  =  11.4  dP> 


out  of  the  detector.  The  TV  product  and  the  signal-to-noise  ratio  into  the 
detector  must  be  sufficient  to  provide  this  SNR  out  of  the  chip  rate 
detector  output  filter  into  the  envelop  detector.! 


5.0  SIMULATION 


5 . 1  Summary  of  Method 

Verification  of  the  analyses  was  performed  through  a  complete 
simulation  of  the  operation  of  the  chip  rate  detector.  Simulated  inputs 
included  single  and  double  fat  differing  rates)  chip  streams  with  Gaussian 
random  noise.  When  two  chip  streams  were  used  as  the  input  they  had 
identical  encoding,  but  were  adjusted  so  that  the  slowest  always  started 
behind  the  higher  rate  stream  to  avoid  corrupting  results  by  the  peaking  at 
coincidence  .  A  block  diagram  of  the  simulated  processing,  using  digital 
signal  processing  techniques,  is  given  in  Figure  5-1. 

Simulation  results  verified  the  analytical  conclusions  within  less  than 
1  dB  over  a  wide  range  of  cases.  Signal-to-noise  ratios  for  single  and 
double  chip  streams,  and  single  and  double  chip  streams  with  noise  found  by 
simulation  agreed  with  the  predicted  ones. 


Input  Filtering 


The  input  filter  to  the  chip  rate  detector  was  varied  over  a  wide  range 
of  bandwidths  before  settling  on  the  method  used  in  this  report.  Simulated 
and  theoretical  results  agreed  for  a  wide  range  of  bandwidths  (greater  than 
the  chip  rate)  as  long  as  the  noise-bandwidth  of  the  filter  was  used  to 
calculate  the  noise  input  power.  A  typical  filtered  waveform  is  shown  in 
Figure  5-2.  In  that  plot,  two  chip  streams  and  noise  are  present.  The  chip 
streams  and  noise  all  have  equal  powers;  namely  0  dB.  signal-to-noise  ratios 
for  both  chip  streams.  As  can  be  seen,  the  chip  streams  are  firmly  embedded 
in  the  noise. 


For  the  purposes  of  the  spectrum  plots  in  this  report,  and  reasonable 
analytical  calculations,  where  a  complete  demonstration  of  the  output 
artifacts  of  the  chip  rate  detector  and  a  verification  of  the  accuracy  of 
the  theoretical  analyses  is  desired,  a  particular  arrangement  of  input 
filtering  was  used.  The  chip  streams  were  unfiltered,  and  the  noise  was 
filtered  with  a  chip  shaped  impulse  response  filter  to  give  a  noise  spectrum 
identical  to  that  of  the  chip  streams.  With  this  choice  of  filtering,  both 
the  predicted  artifacts  (chip  rate  lines  at  higher  harmonics)  and  the  effect 
of  noise  filtering  at  the  chip  rate  bandwidth  can  be  demonstrated.  It 
should  be  emphasized  that  the  basic  results  of  output  signal-to-noise  ratios 
are  not  affected  by  this  choice  (as  verified  by  many  simulation  runs),  only 
the  details  of  the  presentation. 


Filtering 


A  Fourier  transformer  offers  the  opportunity  to  have  both  an  output 
filter  at  the  chip  rate  lines  and  a  spectrum  measuring  device  in  one 
package.  An  integration  time  of  64  chips  (TW  =  64)  was  used  and  a  transform 
length  of  1024.  By  choosing  16  samples  per  chip,  a  chip  rate  of  1/16  of  the 
sampling  rate  was  obtained  with  minimal  aliasing  due  to  fold-over.  Further, 
the  choice  of  the  transform  length  of  1024  was  equivalent  to  a  resolution  of 
1/1024  of  the  chip  rate. 
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To  obtain  smoother  plots  at  t.lio  chip  rate  line,  and  alleviate  the 
sidelobe  leakage  from  transform  channels  at  frequencies  that  did  not  fall 
exactly  upon  one  of  the  transform  frequencies,  cosine-squared  weighting  of 
the  input  data  to  the  transform  was  used  in  the  double  chip  stream  cases. 
With  the  cosine-squared  windowing,  the  effective  TW  product  is  decreased  by 

a  factor  of  2/3rdp  corresponding  to  the  increase  (3/2)  in  noise-bandwidth  of 
the  cosine-squared  filtering.  The  windowing  also  decreased  the  resolution 

by  2/3rds . 

As  the  discrete  transform  of  of  a  real  function  is  redundant  in  its 
second  half,  only  the  first  512  values  (channels)  were  used  in  the  output 
plots. 


5 . 2  Interpretation  of  the  Plots 

For  simplicity  of  interpretation,  the  chip  rate  was  called  100  Hz.  A 
value  of  100  kHz,  or  100  mHz.  could  be  assigned  as  well.  The  Fourier 
transform  interval  between  frequencies  was  1600  Hz./ 1024  =  1.5625  Hz.  With 
the  cosine  squared  windowing,  the  leakage  from  the  discrete  lines  extended 
only  into  adjacent  cells  at  a  meaningful  level.  The  second  chip  stream  was 
chosen  as  110  Hz.  to  ensure  a  good  separation  (resolution)  between  the 
expected  chip  rate  lines.  A  107.  separation  was  well  within  the  analytical 
tolerance  for  verifying  results. 

In  fact,  this  relatively  wide  separation,  compared  to  expected  emitter 
rates,  resulted  in  an  interesting  artifact  that  was  also  predicted  by  the 
analyses.  As  the  chip  rate  detector  has  a  half  chip  delay  only  for  the 
reference  (100  Hz)  chip  stream,  the  110  Hz  stream  does  not  have  a  half  chip 
delay.  The  output  square  wave  for  this  stream  is,  therefore,  asymmetrical 
and  results  in  (small)  even  harmonic  (200  Hz,  400  Hz,  etc.)  chip  rate  lines. 
This  effect  decreases  with  the  square  of  the  percentage  of  asymmetry;  in 
practical  situations  the  artifacts  would  be  undetectable. 

Plots  presented  in  the  following  pages  were  made  by  transferring  the 
output  files  of  the  simulator  (written  in  Turbo  Pascal)  to  Lotus  123 
worksheets.  Use  of  Lotus  123  also  allowed  rapid  and  precise  analysis  of  the 
data.  Signal-to-noise  ratios  of  the  chip  rate  lines  were  calculated  by 
comparing  the  chip  rate  channel  outputs  to  channels  displaced  by  two  where 
only  noise  was  present.  As  the  continuous  spectrum  of  the  noise  is 
relatively  constant  within  this  small  interval,  this  assumption  is  certainly 
val id . 

Additional  checking  of  the  analytical  results  was  performed  to  validate 

2 

the  general  calculations.  Dc-to-chip  rate  line  ratios  (t  / 4  for  a  single 
chip  stream)  were  compared  to  the  theoretical  predictions  and  found  to  agree 

within  less  than  1  dB.  Also,  the  ratio  of  1  to  3  chip  rate  harmonics 
(9)  was  checked;  again  the  results  were  well  within  1  dB. 

To  increase  the  dynamic  range  of  the  accompanying  plots,  continuous 
noise  spectra  were  extracted  from  the  results  of  the  simulation  by  simply 
replacing  the  discrete  line  values  by  the  adjacent  channel  values.  The 
discrete  spectrum  plots  also  include  the  continuous  noise  spectrum.  DC 
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components  were  remove  from  the  discrete  spectrum  plots  by  the  same  method 
to  increase  their  dynamic  range  as  well. 

5 . 3  Statistical  Stability 

Obtaining  good  spectral  estimates,  and  meaningful  signal-to-noise 
ratios  requires  many  runs  of  the  simulation  program.  The  magnitude  squared 
out  of  each  transform  channel  was  averaged  over  400  runs  to  obtain  good 
(stable)  power  spectrum  and  signal-to-noise  ratio  estimates.  At  best,  for 
relatively  independent  transform  channel  outputs,  the  standard  deviation  of 

a  given  channel  would  be  VI / 400  =  1/20^  of  the  expected  value  of  the 
channel  output.  The  two  sigma  value  gives  an  error  of  less  than  ±  107. 
about  907,  of  the  time  (and  greater  the  other  107.  of  the  time).  An 
inspection  of  the  plots  shows  that  the  graphing  adheres  well  to  these 
criteria. 

5.4  Summary  of  Results: 

The  object  of  the  simulation  was  verification  of  the  basic  input-output 
signal-to-noise  equation  (See  Section  2  and  Appendix  VI): 

S2  TV 

0,1  S2  +  S2  +  {a  t2/4)  N2  +  3[S1S2  +/?  SjNj  +  ~0S^\  ' 

with 


S 


Si(2)  =  power  of  chip  stream  1(2), 
and 

Nj  =  noise  power. 

a  and  3  are  chosen,  as  discussed  in  Section  6.1.4,  according  to  the  noise 
spectrum  chosen. 

For  the  chip  shaped  noise  filter,  the  values  of  the  constants  are: 

t2  _  3 
a  1  ~  ?’ 

3=1. 

And  for  the  rectangular  filter  (a  3  pole  Butterworth) ,  the  values  are: 


a 


~1  =  T' 

3  =  m2/ 6. 


o  -  o 


1 


Verification  of  the  equation  was  performed  by  using  numerous 
signal-to-noise  ratios  for  the  input  quantities.  After  a  number  of  trial 
runs  demonstrated  the  correctness  of  the  simulation,  worst  case  runs  were 
made  within  the  limits  of  the  transform  size.  These  runs  are  given  in 
Figures  5-3  through  5-10. 
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Input  values: 


a 


3 

~1  =  7’ 

3  =  1. 


(infinite  input  signal-to-noise  ratio). 


From  the  output  signal-to-noise  ratio  equation  above,  the  theoretical  value 

is : 


rsi 

a 


0,1 


9 


64, 


=  42.67. 

TWa  is  the  product  obtained  after  windowing  (filtering)  of  the  data  prior  to 

transformation  with  the  factor  of  2/3  (due  to  the  noise-bandwidth  increase) 
used  for  the  cosine-squared  weighting. 

From  the  data  which  was  used  to  make  the  plots  of  Figure  5-3,4,  the 
measured  value  was  found  to  be: 


rSlM 


o,l 


42.27, 


which  agrees  remarkably  well  with  the  theoretical  value! 
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Stream  and  Noise,  0  dH.  Sicrnal-to-Noise  Ratios  (Fi 


Input  values: 


S,  =  1, 


S0  =  0, 


Nj  =  2,  (-3  dB.  input  signal-to-noise  ratios) 


2  Q 

^  i r  3 
a  1  =  V 


a  =  i. 


From  the  output  signal-to-noise  ratio  equation  above,  the  theoretical  value 


3  o,l  -  13  a’ 


=  3.28 


From  the  data  which  was  used  to  make  the  plots  of  Figure  5-9,10,  the 
measured  values  for  the  two  chip  streams  were  found  to  be: 


^  Q  1  =  3.22,  (100  Hz.  chip  rate) 


which  is  within  .1  db  of  the  theoretical  value. 
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Double  Chin  Streams.  Infinite  Signal-to— Noise  Ratios 


Inoat  values: 


Sj  =  1, 


S2  =  1,  (0  dB  with  respect  to  reference  chip  stream) 


Nj  =  0,  (infinite  input  signal-to-noise  ratios) 


x2  3 
a  =  7’ 


3  =  1. 


From  the  output  signal-to-noise  ratio  equation  above,  the  theoretical  value 
is: 


SI  _  A  a 
ft  o,l  "  ~o~’ 


=  8.53, 


SI  _  1  a 
ft  o,2  “  T-’ 


=8.53. 


From  the  data  which  was  used  to  make  the  plots  of  Figure  5-5,6,  the 
measured  values  for  the  two  chip  streams  were  found  to  be: 


0  ^  =  8.27,  (100  Hz.  chip  rate). 


ft  0  2  =  (HO  ^z-  rate) 


Both  are  within  .41  db  of  the  theoretical  value,  which  is  quite  good 
considering  the  107.  difference  in  the  chip  rates. 
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Double  Chip  Streams  and  Noise.  0  dB.  Signal-to-Noise  Ratios  fFig.  5-9.10) 
Input  values: 


'I 


“4 


i, 

1,  (0  dB.  with  respect  to  reference  chip  stream) 
1,  (0  dB.  input  signal-to-noise  ratios). 

3 

7’ 

1. 


From  the  output  signal-to-noise  ratio  equation  above,  the  theoretical  value 

is: 


1 

o,l  12.0 
=  3.41 


From  the  data  which  was  used  to  make  the  plots  of  Figure  5-9,10,  the 
measured  values  for  the  two  chip  streams  were  found  to  be: 

M 

1  =  3.36,  (100  Hz.  chip  rate). 

O ,  i 

M 

0  2  =  3.53,  (110  Hz.  chip  rate). 

Both  are  within  .1  db  of  the  theoretical  value,  which  is  quite  good 
considering  the  107.  difference  in  the  chip  rates. 


rs 


rs 
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Single  Chin  Stream  and  Noise.  0  HR.  Signal-to-Noise  Ratio,  Rectangular 


Filter  (Fig.  5-11.12 


Innut  values: 


Sj  =  1, 

S2  =  0, 


Nj  =  1,  (0  dB.  input  signal-to-noise  ratio) 


7T2  X2 

a  ~1  =  T’ 


,5  _  JT 

^  -  T- 

From  the  output  signal-to-noise  ratio  equation  above,  the  theoretical  value 


^  1  TU 

[XJ o , 1  "  ST?  lV 

=  5.08. 

From  the  data  which  was  used  to  make  the  plots  of  Figure  5-11,12,  the 
measured  value  was  found  to  be: 


[Tjo  1  =  4-98,  (100  Hz.  chip  rate). 


The  value  is  within  .1  db  of  the  theoretical  value,  which  is  quite  good. 


ertz 


ertz 


5.4  Conclusions 

Results  of  the  simulation  of  the  chip  rate  detector  confirmed  the 
analyses  and  validated  the  basic  output  signal-to-noise  ratio  formula  given 
in  Section  1  of  this  report  within  1  db. 


6.0  DETAILS  OF  THE  ANALYSIS 


Within  this  section  the  details  of  the  chip  stream  and  noise  modeling 
will  be  described.  Next,  the  output  of  the  chip  rate  detector  will  be 
calculated  for  two  chip  streams  and  noise.  Then,  the  autocorrelation  of  the 
detector,  and  concomitantly,  the  power  spectrum  of  the  the  detector  output 
will  be  found.  The  periodic  chip  rate  components  of  the  output  are 
identified  and  separated  from  the  continuous  (noise)  components.  The  value 
of  the  continuous  spectrum  is  found  at  the  chip  rate  line  as  it  is  the  noise 
component  of  the  signal-to-noise  ratio.  Lastly,  the  signal-to-noise  ratio 
is  calculated  as  the  ratio  of  the  chip  rate  line  (at  a  particular  chip  rate) 
to  the  sum  of  the  total  noise  contribution  from  chip  and  noise  sources  at 
the  chip  rate  line. 

At  best,  the  aforementioned  calculation  is  long  and  tedious. 
Fortunately,  the  result  quoted  in  the  (Section  II)  is  short  and  simple  to 
interpret.  The  basics  of  the  calculation  will  be  presented  in  this  section 
with  repetitive  calculations  of  partial  autocorrelations  (and  spectra) 
relegated  to  the  Appendices. 

6. 1  Modeling 

Definition  of  a  Pseudonoise  Chip  Stream 

A  pseudonoise  chip  stream  x(t)  is  defined  by 

Xj(t)  =  Yj  ck  wl^  "  kTl)’ 
where: 

c^  =  pseudorandom  encoding  sequence, 
w^(t)  =  reference  chip  shape, 

Tj  =  reference  chip  spacing, 

f  _  1 

*1  -  Tp 

=  chip  rate. 

A  rectangular  chip  shape  will  be  used  in  the  final  results,  but  it  is 
unnecessary  to  restrict  the  calculations  at  this  point  to  particulars. 


To  obtain  other  chip  rates  it  is  only  necessary  to  rescale  the  chip 


and  the  chip  stream  at  another  rate  is 


x.2(t)  =  ^  ck  w2(t  -  kT2) , 


=  D  ck  “1 


rT 


1 


-  (t  -  kT2) 
T2 


is  another  chip  stream  at  a  chip  rate  of  f2  =  1  /T2 .  The  factor  of  T^/T2 

also  scales  the  chip  shape,  if  desired,  to  allow  for  the  chip  production  by 
a  clocked  D/A  converter. 

Properties  of  the  PN  Encoding  Sequence 

Statistical  properties  of  the  PN  encoding  sequence  will  be  limited  to 
constant,  unity  amplitude,  zero  mean  sequences  typical  of  spread  spectrum 
systems.  Computations  to  be  encountered  will  require  knowledge  of  the 
second  and  forth  moments: 


Second  Moments: 
E[cjCk]  =  *jk 


A 


1,  j  =  k, 

0 ,  j  ^  k . 


Fourth  Moments: 


E[c.  cc  c  ]  =  6-,  6  +  6-  &  +  6-  &  -  2£-, 

L  k  j  m  nJ  jk  mn  jmu<n  jnaon  jkmn 

with  <5jkmn  =  0  unless  all  the  indices  are  equal  ( j  =  k  =  m  =  n) . 

For  rectangular  chips,  the  fourth  moments  can  be  ignored  as  the  delay 
and  multiply  output  can  be  decomposed  into  constituent  streams  that  avoids 
their  use. 


Multiple  Input  Definition: 

Input  to  the  chip  rate  detector  will  be  assumed  to  be  the  linear 
combination  of  two  different  chip  streams  that  have  the  same  encoding  (but 
different  rates)  and  receiver  noise.  Subsequently,  the  results  for  two  chip 
streams  will  be  found  to  be  easily  generalized  to  many  chip  streams.  And  as 
the  chip  stream  will  be  assumed  to  have  independent  encoding  values  from 
chip-to-chip,  the  results  will  also  hold  for  chip  streams  with  different 
encoding  sequences  (except  when  the  delay  between  the  identically  encoded 
chip  streams  is  less  than  one  chip  time). 
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The  linear  combination  is  given  by: 
s(t)  =  AjXj(t)  +A2X2(t  +  72)  +  n(t) 

where  Xj(t)  and  X2(t)  are  the  chip  streams  defined  above,  is  the  offset 

of  the  second  chip  stream,  and  n(t)  is  the  receiver  noise.  Receiver  noise 
will  be  assumed  to  have  mean  zero  and  an  autocorrelation  function  defined  by 
Rnn(r)  (and  a  spectrum  given  by  the  Fourier  transform  of  R(r)). 

nn 

Output  of  the  Chip  Rate  Detector: 

From  Figure  1,  the  output  of  the  multiplier  is  given  by: 

* 

zm(t)  =  s  (t)s(t  +  Tj) ,  *  =»  complex  conjugation, 

with  rj  =  .5Tp  namely  the  chip  rate  detector  is  designed  to  have  a  half 
chip  delay  upon  the  reference  chip  stream  x^(t).  The  complex  conjugation 
simplifies  later  work. 

In  terms  of  the  constituent  inputs  the  multiplier  output  is  given  by: 

Z„,(t)  =  aJ  X*(t)x,(t  *  rd)  *  x*(t  *  r2)x2(t  *  r2  *  Tj) 

* 

+  n  (t)n(t  +  rd) 

*  * 

+  A!A2  X1 (t)x2^t  +  r2  +  V  +  A1  xi(t)n(t  +  rd^ 

*  * 

+  A2Aj  x2(t  +  r2)x1(t  +  rd)  +  A2  x2(t  +  r2)n(t  +  fj) 

*  * 

+  Aj  n  (t)xt(t  +  rd)  +  A2  n  (t)x2(t  +  r2  +  rd) . 

And  only  four  types  of  cross-product  terms  are  present;  and  this  would  be 
true  regardless  of  the  number  of  chip  streams.  Therefore,  detailed 
consideration  of  these  four  types  of  terms  is  sufficient  to  determine  the 
properties  of  the  expected  output  of  the  chip  rate  detector. 

In  obvious  notation,  these  terms  are: 

I .  Single  Chip  Stream  Output  Product: 

zxx(°’r)  =  Ax  x  (t)x(t  +  r) 
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First,  the  output  chip  stream  product  Type  II  will  be  evaluated  as  Type 
I  is  a  special  case  of  this  product  (x  =  y) .  For  evaluation,  the  chip 
stream  representation  is  inserted  and  the  expectation  over  the  PN  encoding 
sequence  is  taken  to  find  the  average  value.  Performing  these  operations: 


zxy(°,r)  =  AxAy  x  (t)y(t  +  r) 


AxAy  XXCkCj  V*  ”  JTl)  w2(t  ~  kl2  +  7)’ 

k  j 


and  taking  the  expectation: 


E[zxy(°’r>]  *  AxAy  wl(t  ‘  jTl)  U2(t  ■  kT2  *  r)’ 

k  j 


=  AxAy  Xwx(t  -  jTj)  w2(t  -  jT2  +  t) 


If  the  chip  rates  of  the  two  streams  are  identical  (T^  =  T2),  then  the 
expectation  is  a  periodic  function  in  t  with  a  period  of  the  chip  time 


Letting  t  =>  t  +  Tj  and  replacing  the  dummy  summation  index  with  j  =  j  -  1, 

demonstrates  the  periodicity.  Otherwise,  the  function  is  aperiodic  if  one 
chip  rate  is  a  fraction,  close  to  unity,  of  the  other.  Therefore,  the  only 
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chip  rate  lines  are  the  individual  lines  that  would  come  out  of  the  chip 
rate  detector  had  only  single  chip  streams  been  present. 

The  cross  product  terms  will  only  contribute  when  they  line  up  at  some 
time  given  by  solving  the  simultaneous  equations  (found  by  equating  the 
arguments  of  a  particular  chip  to  zero): 

t  -  jTj  =  0, 
t  -  jT2  +  r  =  0 
or, 

-  T2)  =  -r, 


At  this  chip  time  the  two  different  chip  streams  will  overlap  and  produce  a 
"blip"  that  lasts  for  a  time  depending  on  the  difference  in  their  rates. 
"Blip"  time  can  be  found  by  similar  equations  at  the  leading  and  trailing 
edge?  of  chip  line  up.  According  to  their  range  (r),  and  chip  rates,  this 
coincidence  can  have  occurred  in  the  past  or  will  occur  in  the  future. 


By  simple  examination  of  the  expectation  of  the  expression  above,  for 
equal  chip  times,  the  frequency  content  of  the  discrete  lines  may  be  found. 
A  general  Fourier  transform  is  not  necessary  to  see  the  result  for 
rectangular  chips.  Consider  the  expression  above  for  x  =  y: 

E[Zxx(°,r)]  =  A^  £wl(t  -  jTj)  Wj(t  -  jTj  +  t)  . 
j 

With  t  -  Tj/2,  a  half-chip  overlap,  the  summation  can  be  seen  ,  graphically, 

to  be  equivalent  to  a  simple  square  wave  at  the  chip  rate  as  in  the 
following  Figure  6-1.  Of  course,  this  was  the  reason  for  use  of  a  chip  rate 
detector  in  the  first  place. 


Figure  6-1  Chip  rate  Detector  Output 
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Chip  rate  lines  will  appear  at  odd  harmonies  of  the  ehip  rate  plus  a  DC 

term  of  .5A^  as  expected  of  chip  rate  detectors.  A  transform  of  the  general 

expression  above  can  be  performed  to  provide  the  line  amplitude  for 
arbitrary  chip  shapes;  but  this  operation  will  be  deferred  until  after  the 
autocorrelation  function  of  the  output  is  calculated  to  find  both  the 
discrete  and  continuous  spectra. 

For  a  chip  stream  other  than  the  reference  stream,  the  square  wave  is 
asymmetrical  due  to  the  fact  that  the  refeience  half-chip  delay  is  not  half 
a  chip  in  this  instance.  The  spectrum  will  then  include  both  odd  and  even 
chip  rate  harmonics  terms  at  the  chip  rate  of  the  particular  stream  under 
discussion. 

The  increasing  separation  of  higher  harmonic  lines  from  the  reference 
chip  stream  may  be  a  possible  identifier  at  a  sufficient  signal-to-noise 
ratio. 


As  the  noise  is  independent  of  the  chip  stream,  and  has  zero  mean,  this 
type  of  term  also  has  a  zero  mean,  and  therefore,  does  not  contribute  to  the 
mean  output. 


The  expectation  of  this  term  is  straight  forward  and  easy  to  interpret. 
Taking  the  expectation  of  znn(t»r): 


E[znn(°»7')]  =  E[n  (t)n(t  +  r)] , 


for  stationary  noise.  The  expected  mean  value  of  the  noise  is  a  DC  term 
which  is  merely  the  input  noise  autocorrelation  evaluated  at  one-half  chip 
time  offset  from  the  peak.  And  this  is  to  be  expected;  a  zero  delay  (r  =  0) 
would  simply  be  a  power  measurement  of  the  noise. 


Summary  of  the  Chip  Rate  Detector  Output  Expectations 


Expected  outputs  of  the  chip  rate  detector  are  what  might  be  expected. 
Discrete  spectral  lines  of  each  chip  stream  are  present  and  depend  upon  the 
chip  shape.  For  a  rectangular  chip  stream,  the  reference  chip  stream  has  a 
symmetrical  square  wave  as  its  output;  precisely  what  would  be  obtained  if 
it  were  the  only  input.  Other  chip  streams,  with  differing  chip  rates,  have 
asymmetrical  square  waves  (as  the  half-chip  delay  is  not  half  a  chip  for 
rates  differing  from  the  reference  chip  stream)  at  their  own  chip  rates.  As 
these  waves  are  asymmetrical,  they  also  have  even  harmonics  present. 
However,  these  will  be  quite  small  for  chip  rates  close  to  that  of  the 
reference. 
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Mixing  of  chip  streams  produces  a  "blip"  in  the  output,  when  the  chip 
streams  overlap;  this  output  can  be  in  the  present,  past,  or  future 
depending  upon  range  and  chip  rates  and  start  of  the  observation  time. 

Mixing  of  chips  and  noise  produces  a  mean  zero  output  with  no  periodic 
content . 

And  as  expected,  the  noise  output  expectation,  a  DC  voltage,  is  no  more 
than  the  noise  autocorrelation  function  evaluated  at  a  half-chip  time. 

6.3  Calculation  of  the  Chip  Rate  Detector  Output  Spectra 

Expected  values  of  the  detector  output  have  given  insight  into  the 
operation  of  the  chip  rate  detector.  The  output  spectrum  allows  the  output 
signal-to-noise  ratios  to  be  developed.  Evaluation  of  the  output  spectrum 
of  the  chip  rate  detector  is  straightforward,  but  arduous.  The  simplest 
approach  is  the  calculation  of  the  time  averaged  autocorrelation  as  the 
detector  chip  stream  outputs  are  not  stationary.  This  generalized  quantity 
also  allows  a  generalized  spectrum  to  be  calculated  as  in  the  stationary 
case  (which  is  a  specialization).  They  are,  however,  cyclostationary  which 
allows  the  generalized  spectrum  to  identify  both  the  discrete  and  continuous 
spectral  components. 

Calculation  of  the  Chip  Rate  Detector  Autocorrelation  Function 

First  the  output  of  the  delay  and  multiply  will  be  evaluated.  The  time 
averaged  autocorrelation  function  is  defined  by: 


R"z(r>  =  ‘im  CT  f  E[Z*(t)2B(t  ♦!■)]  dt. 

1  “*00  J  ^ 


For  future  use  the  time  averaged  autocorrelation  function  will  be  hereafter 
simply  called  the  "autocorrelation  function;  and  the  operation  will  be 
identified  by: 


T 

E[*]\  =  lim  ^  f  E [ - ]  dt. 

'  T-ao  J  T 


And  using  the  multiplier  output  from  above: 

Rzz(f)  =  (E[s(t)s  (t  +  rd)s  (t  +  r)s(t  +  rd  +  r) 
with  s(t)  given  again  bv: 

s(t)  =  Aj  Xj(t)  +A2  X2(t  +  +  n(t). 


**.<*.<  of oc  on 


Tn  the  previous  notation,  where  the  four  types  of  cross  products  were 
identified,  the  delay  and  multiply  output  can  be  written  a  little  more 
simply: 

zm(t)  =  s  (t)s(t  +  rd)  j 

I 

l 

=  zn(0,rd)  *  z22(r2,r2  *  rrf)  ♦  znn(0,rd)  (a) 

*  z,2(0,r2  *  rd)  *  z21(r2,rd)  (b) 

*  zln(0'rd)  *  z»l(#’rd>  (c) 

*  z2n<’2’rd>  +  zn2^0,r2  +  rd)-  <d> 

Next,  the  correlation  products  must  be  formed,  which  is  tedious  at 
best.  To  reduce  the  number  of  terms,  a  simplifying  assumption  will  be  made. 
It  will  be  assumed  that  the  two  input  streams  never  line  up  during  the 
period  of  observation;  the  "blips"  identified  in  the  multiplier  output  above 
are  not  present.  In  this  case,  the  codes  may  be  considered  as  independent. 
In  the  analysis,  the  code  X2(t)  will  have  a  code  sequence  b^,  rather  than 

Cj,.  otherwise,  the  code  properties  will  remain  the  same. 


6.3.1 


;ntiti( 


iut  Autoc 


)n  lerms 


Fortunately,  by  virtue  of  taking  expectations,  many  of  the 
cross-product  correlation  terms  will  vanish.  For  example,  the  expectation 
of  the  product 

E[ZH  (°,  rd)zi2(0,  Td  +  = 


k  j  m  n 


•Wj(t  ~  jTj)  Wj(t  -  kTj  +  rd) 


1  w< ( t  -  jT.  +  r)  w9(t  -  kT9  +  r.  +  r) , 


is  zero  as 


E[ckcjc«bJ  =  EtckcjcJ  Etbn3 


for  the  independent,  mean  zero  coding.  Any  correlation  product  containing 
only  single  members  of  the  coding  streams,  or  noise,  must  vanish  in  a 
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similar  manner.  No  cross-products  between  members  of  lines  (a)  through  (d) 
will  be  retained! 


The  remaining  terms  are^ 

zm(t)  =  s  (t)s(t  +  rd) 

* 

=  zn(°’rd)zll(r,rd  +  T>) 


*  z22<r2-r2  *  rd>z22<r2  *  r'r2  *  rd  *  r> 

*  z^(°’rd)zim(r’rd  +  r> 


*  2ll(0,rd)z22(r2  *  r,r2  *  rd  +  r) 

*  z22 ( r2 ' r2  *  rd>zll(r’rd  +  r> 

*  zll^0,7"d>znn^r’  rd  *  r) 

*  znn(°'rd)zll(r’rd  *  r) 

*  z22<r2’r2  *  rd>  z„i,(r'rd  +  r> 

*  znn(0'rd)z22<T2  *  r’r2  *  rd  *  r> 

*  z12^°’r2  *  rd>z12<r’r2  *  rd  *  r> 

*  z12^0, r2  *  rd)z21<r2  *  r’rd  *  r) 

*  z21^r2,rd^z12^r,r2  *  rd  +  r> 

*  4<r2’rd>z21<r2  *  r’rd  *  r> 


*  zli.<°'rd'zln<r’rd  f  r> 


I 


zl,.(0'rd)znl(r'rd  '  r|' 


*  znll°'rd)zln<r'rd  *  r) 


+  zr,l(°'rd>znl<r’rd  *  r> 


*  z2n<r2’rd>  z2„<r2  *  r'rd  *  r> 

*  z2n<r2’rd>zn2<r-r2  *  rd  +  r> 

*  zn2(0’r2  *  rd)z2n(r2  *  r'Td  *  r) 

*  zn2<0-r2  *  Td^zn2^T,r2  +  rd  +  r~> 

Notation  for  the  individual  terms  may  again  be  simplified  by  writing: 

zxyuv‘rl,r2’r3’r4>  =  zxy<rl’r2>  zuv<r3-r4> 

A  *  * 

=  AxAvAuAv  +  rl)y  ^  +  r2^u  ^  +  r3^v^  +  tA) 


with  the  partial  autocorrelations  of  each  term  denoted  by  a  similar  notation 
(1  «  x,  2  «  y,  n  m  n): 

Rxyuv<rl,r2'r3'r4>  -  (E  [zxy<rl'r2>  zuv<r3'r4>])- 

Using  this  notation,  the  delay  and  multiply  output  autocorrelation  is: 


Oz>  =  (E  zm<1  *  r>l) 

=  Rmi<°'rd’r'rd +  T) 


(Individual  Input  Terms) 


+  R9999(r9  ’ Tc)  +  rH’r9  +  T'  r2  +  rd  +  r) 


2222 v  2  ’2  d’  2 
*  Rnnnr,<°'rd'r'rd  * 


*  R1122(0’rd’r2  *  r,r2  *  ;J  *  r) 
+  »2211(r2’r2  +  rd’z’rd  *  r) 


(DC  Terms) 
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+  Rll1m<0'rd’r'rd  *  r) 

*  RMll(°'rd'r'rd  *  r) 

*  R22nn(r2’r2  *  rd'r'rd  *  r> 


*  Rnn22<°’rd'r2  *  T’Tr>  *  r> 


2  d 


*  R1212<°’r2  *  rd'r'  r2  *  rd  *  r) 

*  R1221<°’r2  *  rd'r2  *  T'Ti  *  T) 

*  «2112<r2’rd'r'r2  *■  rd  +  r> 

*  R2121(r2’rd’r2  *  r’rd  *  r> 


(Chip  Mixing  Terms) 


*  Rlnln<°'rd'r'rd  *  r> 


(Noise/Chip  1  Mixing  Terms) 


*  Rliml<°’rd’r’rd  *  r> 


*  Rnlln^°’rd’r,rd  *  r> 

*  Ri,lnl(°'rd>T’rd  *  r> 

*  R2im2(r2>rd'r’rd  *  r> 

*  R2nn2(r2’rd’r'r2  *  rd  +  r> 

*  Rn22n<°'r2  +  rd’r2  *  r’Td  *  r> 

*  Rn2n2^0,rd,r,r2  *  rd  *  r) 


(Noise/Chip  2  Mixing  Terms) 


6.3. 2  Evaluation  of  the  Individual  Autocorrelation  Terms 


All  the  individual  terms  do  not  have  to  be  evaluated  separately;  only 
the  few  classes  identified  above.  Individual  input  terms,  those  that  would 
be  present  if  no  other  input  be  present,  will  be  evaluated  in  detail  as 
exemplar  cases.  Details  of  other  terms  will  be  relegated  to  Appendices. 


Individual  Input  Terms 


These  terms  are  present  with  only  a  single  input.  Other  types  of  terms 
deal  with  the  interaction  between  terms. 


It  is  sufficient  to  compute  the  autocorrelation  product  R9900(r); 
is  a  special  case  with  r9  set  to  zero.  In  detail,  the  product  is: 

R2222(r2,72  +  rd'r2  +  r'  r2  +  rd  + 

=  <a4  I  III  E[ckcjV„] 

k  j  m  n 

•w^t  "  jT-2  +  r2)w2t  -  kT2  +  +  r^) 

•w9(t  "  J^2  +  r2  +  r)w2^  "  ^2  +  r2  +  rd  + 

Computation  of  this  autocorrelation  can  be  performed  analytically  using  the 
fourth  moments  of  the  encoding  sequence  for  an  arbitrary  chip  shape;  but  it 
is  simpler,  and  more  intuitive,  to  use  graphical  methods  with  rectangular 


Figure  6-‘2b.  Deterministic  Product  Chip  Rate  Waveform 


Figure  6-2c.  Random  Product  Chip  Waveform 


From  Figure  6-2  the  output  of  the  delay  and  multiply  can  be  seen  to  be 
composed  of  two  waveforms.  A  deterministic  square  wave,  with  a  period  of 
the  chip  rate,  occurs  due  to  the  same  chip  being  multiplied  by  itself  for  a 
time  rQv  =  as  in  Figure  6-2b.  A  random  waveform,  as  in  Figure 

6-2c,  develops  in  the  region  between  overlaps.  The  random  waveform  has 

equilikely  ±A  states  which  have  a  duration  of  Tj.  Also,  the  random 

waveform  is  completely  uncorrelated  with  the  deterministic  waveform  as  each 
is  zero  where  the  other  has  a  value. 

Therefore,  the  autocorrelation  of  the  composite  waveform  is  the  sum  of 
the  autocorrelations  of  the  waveforms  in  Figure  6— 2b  and  6-2c,  and  is  given 

by : 


4  r 

^xxxx^  "  ^  L^(Unit  Square  Wave  Autocorrelation)^ 

+  ^(Random  Biphase  (Duration  rj,  Interval  T9))^r^’ 


I 

m 


t  tut*  « re-r^TT^n  -v  r^r  ^n 


whore  from  Appendix  T: 


^(Unit  Square  Wave  Autocorrelation)^ 


*  *T  1 

k 


c \  ikii,’0r 

e  z 


with  the  square  wave  coefficients 


1/2,  k  =  0, 


fk  =  (_1)  ik’  k  oddl 


0,  otherwise. 


i.c  =  o,’c  =  2t/Tc  is  the  chip  rate. 


As  the  spectrum  is  the  transform  of  the  autocorrelation  function,  it 
has  only  discrete  lines  at  DC  and  the  harmonics  of  the  chip  rate.  The 
spectrum  can  be  written  in  terms  of  impulse  functions: 


SSVM  =  2 7r  A4  £f£  6(lj  -  ku,-c) 
k 


2 

=  A4  <5(0)  +  ^  S(u  -  ku;c) ,  1  =*  exclude  k  =  0. 


C  r 

The  value  of  the  1  chip  rate  line,  the  signal  output  component  of  the  chip 
rate  detector  is,  therefore: 


And  from  Appendix  II, 


"(Random  Biphase  (Duration  r^,  Interval  T2))^r^’ 


=  A4  y  J  Wd(t)wd(t  +  T)  dt, 


C  — oo 
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where  w ^ ( t )  is  the  rectangular  chip  of  width  (Figure  6-'2c),  and  the 

autocorrelation  is  a  simple  triangle  with  a  base  of  2r^  and  a  height  of 
-j/T.-).  As  the  spectrum  is  the  Fourier  transform  of  the  autocorrelation 
'Appendix  II): 


,  t  4  4  sin (ur^/2) 

=  *  TJ  -rd/2  • 

,  T0  sin(u/T0/4)  2 

=  4  — 7T— n -  »  for  T,  2  To/2 , 


which  is  a  power  spectrum  that  is  twice  as  wide  as  the  input  chip  stream 

spectrum. 

At  harmonics  of  the  chip  rate  where  a  narrowband  filter  will  be  placed 
to  detect  the  chip  rate  line  the  continuous  spectrum  will  have  the  value: 


T 

-dr ' "']< )  =  ^  ~ — 7  sin2(mk/2) 

W 


( b . )  Noise 

Development  of  the  delayed  and  multiplied  noise  autocorrelation 
requires  assumptions  about  the  statistics,  and  spectrum,  of  the  noise.  The 
bandwidth  of  the  noise  and  it's  spectrum  do  not  change  the  basic  physics; 
but  they  do  change  the  exact  form  of  the  output  signal-to-noise  ratio 
constant  and  the  difficulty  of  the  calculations.  The  autocorrelation  of  the 
multiplier  output  was  found  to  be: 

Rnnnn(r)  =  (Etn(t)n  ^  +  r<pn  ^  +  r)n(l  +  rd  +  r^)’ 

=  E[n(0)n  (rd)n  (r)n(rd  +  r)] , 

=  EPWJ 

assuming  the  noise  process  is  stationary. 

From  J.  (1.  Laning,  Jr.  and  R.  II.  Battin,  "Random  Processes  in  Automatic 
Control",  McGraw-Hill  Book  Company,  Inc.  1956,  the  expectation  can  be  found 
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«*_  s .  s  _  ■'  j' 


for  Gaussian  processes  as 


E[X1X2X3X4]  -  m12m34  +  ml3m24  *  m14m23. 


where 


'i j  =  E[x(r;)x(rj)] 
=  R*x<ri  *  rj)' 


In  this  instance  the  correlation  R  (r)  is  given  bv: 

nnnnv  '  ° 

R  (r)  =  R  (r.)2  +  R_  (r)2  +  R__(r  +  r.)R  (r  -  r.), 
nnnnv  '  nnv  d'  nnv  '  nnv  d'  nnv  d; 

where  R  (r)  is  the  autocorrelation  function  of  the  input  noise  process. 

The  output  noise  spectrum  is  obtained  through  taking  the  transform  of 
Rnimn<r>: 


S  (u) 

n  n  rtn  V  • 


=  2?r  Rnn(rd)2  +  h  /da;l  Sm)(a;l)SnnK  '  ^ 


i  C  i2u<  ta 

+  I5F  J  da;l  Snn(tJl)Snn(a;l  “  e 


where  Snn(u;)  is  the  spectrum  of  the  input  noise  given  by: 


S»  =  f  R  (r)  e~luT  dr. 
nn'  '  J  nn'  ' 


For  computational  purposes,  two  autocorrelation  functions  will  be 
chosen  for  the  noise,  or  equivalently,  spectra. 

I •  Noise  Using  Chin  Spectra 

Noise  spectra  may  be  chosen  to  be  identical  to  the  chip  spectra: 
namely,  having  the  same  basic  bandwidth.  In  this  case,  the  noise 
autc :orrelation  function  is  given  by  the  usual  triangle: 


R  (r)  =  Nt 
nnv  '  I 


[1  -  .  [t |  <  Tj , 

0  ,  otherwise. 
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*  V 

M  w  m 


"VVVV'! 


with  a  spectrum 


snn(^)  =  Ni  Ti 


sin(o;T1/2) 


As  this  is  identical  to  the  spectrum  of  a  chip  stream,  the  formula  of 
Appendix  IV  for  and  S^n(u;)  may  be  used  for  the  spectrum 

Sr*).  Therefore, 
nnnnv  ' 


=  /  Ni2  m 


*  Ti  h 


2  2 
1  (^T,); 


[2  -  cos(a)T^/2)]  1  +  2 


sin(u;T1/2) 


At  u  =  a;,  =  2^/T,,  the  chip  rate,  the  multiplier  output  spectrum  amplitude 
due  to  the  noise  component  is: 


W“i>  ■  1  T1  ‘V  h 


II .  Rectangular  Noise  Spectrum 

A  rectangular  noise  spectrum,  of  a  width  of  twice  the  chip  rate,  may 
be  chosen.  It,  however,  complicates  later  chip/noise  calculations.  The 
rectangular  spectrum  is  given  by: 


’-7-  Nr,  H  <  w  , 

s„„(-)  =  1  c  , 

[  0,  otherwise. 

From  earlier  work,  the  multiplier  output  spectrum  is  given  by: 

Snnnn(u:)  =  2t  6(0)  Rnn(r,)2 
nnnnv  '  v  '  nnv  d' 

1  r  Tl2  „2  .  >  f,  .  sil>[(2u,l-“')r(1]1 

Ti  NI  <2jl  1  +  (2^-u/J  rd - 


*  ii  Ni  <2“i  [>  * 

for  <  2u^,  and  is  zero  otherwise. 
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At  the  chip  rate  the  multiplier  continuous  noise  spectrum  output  has 
the  value: 


W“'1>  ■  T1 1  Nr 


This  value  will  produce  the  classical  results  for  the  output  of  a  chip  rate 
detector  for  a  single  chip  stream  and  noise.  But,  this  spectrum  makes  the 
rest  of  the  work  in  this  paper  extremely  difficult.  The  triangular  noise 
autocorrelation  above,  with  the  bell-shaped  spectrum,  allows  the  rest  of  the 
analysis  to  be  much  simpler. 

And  the  difference  in  the  output  noise  powers,  of  the  two  cases,  is 
easily  calculable.  Taking  the  ratio  of  the  cases  I  and  II: 


<tT7’ 

=  .608, 

=  -2.12  dB, 

Hereafter,  the  noise  autocorrelation  of  case  I,  the  triangular  form,  will  be 
used  in  the  analytical  calculations  with  the  results  verified  by  simulation. 

In  the  final  results  the  noise  contribution  constant  will  be  identified 


C  =  a 


Then,  for  comparison  with  other  chip  rate  detector  work,  a  can  be  chosen  as 
unity  or  the  value  given  above  for  simplified  formulae  (Cn  =  3/2). 


From  the  efforts  of  the  previous  section  and  the  Appendices,  the  total 
discrete  and  continuous  contributions  to  the  chip  rate  detector  output  can 
be  collected.  Work  in  this  section  has  identified  four  types  of  output 
terms  from  the  chip  rate  detector.  These  types  and  the  pages  upon  which 
they  may  be  found  within  the  derivations  are: 

I .  Chip  Rate  Sine  Waves: 

An  individual  chip  rate  line  power: 


>  =-4 
1  2’ 


(Page  6  -  14) 


9  t  h 

Sj  =  Aj,  the  input  power  of  the  j  n  chip  stream  (at  baseband). 

II.  Chip  and  Noise  Independent  "Noise"  Spectra  Outputs 
Individual  chip  stream  continuous  spectrum  value  at  the  chip  ra 


>C  ,2  4 

j  j 


(Page  6  -  15) 


1/Tj  =  chip  rate  of  the  j  n  chip  stream. 

Noise  continuous  spectrum  value  at  the  chip  rate: 
Chip  shaped  impulse  response  filter: 


’nnnn^  =  1  T1  NI 


(Page  6  -  17) 


Rectangular  noise  spectrum: 


,H>  ■  h  1  "i, 
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III.  Chip/Chio  and  Chip/Noise  Intermodulation  Spectra 
Chip/chip  mixing  continuous  spectrum  value  at  the  chip  rate: 


Pjk  =  T!  Sj  Sk  4’ 

j  j  7T 


(Appendix  -  14) 


with 


1/Tj  =  chip  width  of  the  chip  streams  all  considered  approximately 
equal . 

Chin/Noise  continuous  spectrum  value  at  the  chin  rate: 

Chip  shaped  impulse  response  filter: 


Snnnn(a,l)  =  T1  SjNI  T 

«  IT 


(Appendix  -  18) 


Rectangular  noise  spectrum: 


S  (uO  S  N. . 

nnnnv  V  2  1  j  1 


(Appendix  -  19) 


6.3.4  Calculation  of  the  Chip  Rate  Detector  Output  Signal-to-Noise  Ratio 

By  combining  the  various  terms,  above,  the  output  signal-to-noise  ratio 
may  be  found. 

First,  the  chip  rate  sine  wave  power  for  the  j*'*1  chip  stream  (I)  is 
identified  as  the  signal  power  output.  It  is 


■4 


Then,  the  continuous  chip  stream  spectrum  powers  (II  and  III),  at  the  chip 
rate  are  summed  together,  multiplied  by  the  narrow  band  chip  aetector  output 
filter  bandwidth  (B  =  1/T),  to  form  the  total  noise  output  (for  two  chip 
streams) : 


No  =  B  “2  [  s?  *  sl  *  7  NI  *  3  tSl  S2  *  SjNj  ♦  S/j  ] , 
and  the  chip  detector  output  output  signal-to-noise  ratio  is: 


NjoJ  "  FT’ 


-1  .  '  T  -  '  5T  -  — — 

T  >  1  '  '■  &T‘-  TV 
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■  ,ij,miWiimBi„u,WB  niuuamnaiw. 


TW 


3,2., 

7  NJ  3 


[S,  S2  ♦  S,Nj  *  S2Nj 


with  V  =  1/T^  as  the  chip  rate  in  hertz,  and  T,  the  bandwidth  (or 
integration  time  of)  the  output  filter. 


Rectangular  Noise  Spectrum 


No  = 


T1 


C2  .  C2  .  7T 

bl  b2  T 


N?  *  3 


[SlS5 


[S,NI  *  S2Nj] 


and  the  chip  detector  output  signal-to-noise  ratio  is: 


TW 


S?  +  S2  +  T  N?  +  3  [S1S2  +  [SjNj  +  S2Nj] 


,  » 


And  we  have  the  maior  result  of  this  report! 

By  defining  two  constants  a  and  /?,  the  two  cases  (chip  shape  and 
rectangular  input  noise  spectra  can  be  written  in  a  single  form  for  use  in 
the  rest  of  this  report.  Define: 

Chip  Shaped  Input  Noise  Spectrum: 


and 

a  =  1; 

Rectangular  Input  Noise  Spectrum: 

Or  =  1, 


With  these  values,  the  signal-to-noise  ratio  for  both  cases  both  cases  may 
be  written  in  the  form: 


V1 

SjS 
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[  +  a  \  NJ  +  3  [SjS2  +  0  [SjNj  +  SjNj] 


which  is  the  form  used  throughout  this  report. 
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APPENDIX  I 

Unit  Square  Wave  Autocorrelation 


Expansion  of  the  periodic  chip  function  in  a  Fourier  series  is 
relatively  simple.  An  expansion  is  given  by: 


5 

l'. 

c~t 

II 

yi 

i 

fk  e 

V" 

*  , 

with  coefficients 

I 

JJ2 

p 

K 

£ 

II 

M-h 

f 

-V2 

— i  2  Trkrp — 

f(t)  e  Lc  dt. 


For  the  unit  square  wave  of  a  width  of  half  a  period,  the  coefficients  are: 


1/2,  k  =  0, 

.k  1  t.  ,ji 


fk  =  -<-»*•  k°d(i' 

^  0,  otherwise. 

The  autocorrelation  function  is 


•here  =  2mk/Tc  is  a  harmonic  of  the  chip  rate. 


As  the  spectrum  is  the  transform  of  the  autocorrelation  function,  it 
has  only  discrete  lines  at  DC  and  the  harmonics  of  the  chip  rate. 


sss 


A  unit  amplitude,  binary,  random  square  wave  can  be  represented  as 


f(t)  =  X  Ck  "  kTc)- 

The  time  averaged  autocorrelation  function  is  given  by 

,T 

Rff(r)  =  lim  ^  J  E[f(t)f(t  +  r)]  dt, 

T~*oo  -T 

J 


=  lim  7T  /T  X  E[ckcj^  "  VcW  ~  jTc  +  r)  dt’ 
T"°°  "T  j  i  k 

rT 

=  lim  J  V  w(t  -  kTc)w(t  -  kTc  +  r)  dt, 

T-*oo  -T 


as  E [c^c j]  =  for  uncorrelated  binary  sequences  with  values  of  ±1. 

To  evaluate  the  expression,  change  the  variable  of  integration  and  the 
order  of  summation.  The  resulting  expression  is 

T-kT 


Rrr(r)  =  lim  1  Y  [  w(t)w(t  +  r)  dt, 

11  T-vr,  4*  •'-T-kT 


T'*00  k  -T-kTc 


00 

=  y  j  w(t)w(t  +  r)  dt, 

C  —oo 

upon  letting  T  =  NTC  and  taking  the  limit  on  N,  and  assuming  the  chip  shape 

falls  of  reasonably  rapidly  (which  it,  of  course,  does  for  rectangular 
chips). 

For  rectangular  chips  of  width  T  ,  the  autocorrelation  is  a  simple 

triangle  with  a  base  of  2Tw  and  a  height  of  Tw/Tc-  As  the  spectrum  is  the 

Fourier  transform  of  the  autocorrelation,  the  simple  relation  for  the 
spectrum  is: 


using  simple  transform  relations,  and 


— OC 


is  the  Fourier  transform  of  the  chip  shape. 

For  the  rectangular  chip  of  width  T  ,,  the  power  spectrum  is 
w  w 


w 

r 


sin(<jTw/2) 


2 


1 


which  is  the  typical  spectrum  for  rectangular  chips.  The  spectral  width  to 
the  first  null  is  the  radian  frequency  =  27t/Tw .  Therefore,  a  chip 

rate  detector  which  has  a  delay  of  one-half  chip  has  a  continuous  power 
spectral  output  component  that  is  twice  as  wide  as  the  input  chip  stream  )as 
might  be  expected) . 

At  the  chip  rate,  the  self-noise  power  spectrum  out  of  the  chip  rate 
detector,  of  a  single  chip  stream,  has  the  value  (assume  Ty  =  .5  T£) : 
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Expansion  of  Periodic  Summations 


Expansion  of  the  periodic  chip  summations  in  a  Fourier  series  is 
relatively  simple.  As  before,  an  expansion  is  given  by: 


_  127TKJ- 

f(t)  =  £  fk  e  c  , 


with  coefficients 


=  T —  J 

k  *c  -Tc/2 


Lc/Z  — i  2  7rkrp- 

f(t)  e  V  dt. 


and  where 


f(t)  =  ^  w(t  -  kTc)w(t  -  kTc  +  r). 
k 

Inserting  the  summation  in  the  coefficient  expression: 


rlc'2 

'k  =  T~  J  ^  w(t  -  kTc)w(t  -  kTc  +  t)  e  ‘c  dt. 
c  -Tc/2  k 


-i2xkn 


,Tc/2  -  «c 


.  r'C~  ”'c  -i27rkm:- 

=  y-  >  I  w(t)w(t  +  r)  e  c  dt. 

Tc  f  ~Tc/2  -  kTc 


And  finally: 


00 

fk  =  j-  /  w(t)w(t  +  r)  e 


-i2^kn 


c  “oo 
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DC  k  Chip  Mixing  Terms: 

These  terms  only  occur  due  to  the  interaction  of  chip  streams  with  each 
other,  or  with  noise,  entering  the  chip  rate  detector. 

Chip  Streams 

A  typical  chip  mixing  term  is 

Rll22(rr  r2,r3,r4^  =  (  A1A2  E^CkcjlEtbmbJ 

k  j  m  n 

* 

" w t ( t  -  jT1  +  r{)  wt(t  -  kT 1  +  r2) 

•w.2(t  -  jT2  +  7*3)  w2^  _  kT'2  +  r4^) 

=  (  A1A2  X  Wl(t  '  jTl  +  rl}  Wl(t  ~  jTl  +  r2} 

k 

•  £  w*(t  -  kT2  +  r3)  w2(t  -  kT2  +  r4)) 

j 

Other  chip  mixing  terms  are  merely  changes  of  subscripts  and  times  (r  1 s); 

they  are  different  cases  of  the  correlation  product  under  consideration. 


Now,  the  summations  are  periodic  functions  in  the  chip  times  and  can  be 
expanded  in  Fourier  series.  For  example: 


A 


an  e 


inu^t 


*1  = 


fi(t>  -  I 

=  £  Wj(t  -  jT1  +  T{)  w*(t  -  jT1  +  r2) 


and  from  Appendix  III: 


a  = 


ina;.  r.  .  *  -inaut 

e  j-  J  wj(t)  wj(t  +  r2  “  ri)  e  1  dt. 

1  — 30 


As  the  chip  shape1  is  given  in  forms  of  its  transform  by 


00 

=  h  /  Wl^’ 


\  i  u/t  i 

)  e  du 


m^,  r 


an=  e 


ri 


1 

Tz 


^2 


.00  *00  ^00 


T7  /  dt  J  du/'J  dw'  V^u/'jV* 


1  —oo  —oo  -00 

.ei(ui"-u'-nu1)t  e-iu-(  r2  -  r,)  dt_ 


(V  > 


in  a;.  r. 


.00 


a_=  e  1  1  T"  /  dw'  Wl(u;'  +  ^  e~itJ' (  ^  ”  ri)- 

1  -oo 


n 


Similarly,  the  other  periodic  function  can  be  written  as  a  Fourier 

series : 


f  m  A  V  A  eimw2t  ,  _ 
f  2  ( t )  —  cl__  e  •,  u)ci  — 


m 


2 


=  £  w2(t  -  kT2  +  r3)  w*(t  -  jT2  +  r4) 
k 


with 


.00 


V  dw"  *  -“PV*'”)  e~iw"(  T4  - r3>- 


-oo 


Finally,  the  time  averaged  expectation,  using  the  Fourier  series 
representations,  is 


9  9  V1  *  /  Hna,i  imu;2)tv 
R1122(rl,r2,7"3,r4^  =  A1A2  andm  \e  / 


m  n 


Time  averages  of  sine  waves  vanish.  As  w^/u/2  is  close  to  unity  (clock 

rates  are  similar),  the  only  contributor  of  interest  is  the  zeroth.  Higher 
order  contributors  will  occur  at  Nu^  =  Mu;2  if  -  M/N.  When  M/N  is 

close  to  one,  these  contributors  occur  at  frequencies  where  the  chip  rate 
lines  are  too  far  down  to  be  of  interest.  Therefore, 


Rll22^rl,r2’r3’r4^ 


y ill** ly 2y iv iy .-ov  v.'~ .-~>yv  .-'■••y- 


6 


.2*2  0  ,* 

=  A1A2  aOdO 


=  A!»2Rr/  <*"'  |W,(W')|2e-iu'(r2-rl). 

1  — 00 

•  ^  J  dw"  |V*(u/')|2  e  iw"(r4  "  r3) . 

2  — 00 

Other  pair  products  follow  an  identical  development.  They  are  simply 
-earrangements  of  the  arguments  r  .  For  example,  the  term 

l12l2^rrr2’r3’r4^  =  (  A1A2  XXSX  E^ckcnJE[EjbJ 

k  j  m  n 

•Wj(t  -  kT1  +  rj)  w2(t  -  jT2  +  r2) 

•w*(t  -  jT1  +  r3)  w2(t  -  kT2  +  r4)^ 

=  ^  42A2  wt(t  -  kTj  +  rt)  w2(t  -  jT2  +  r2) 
k 

•  ^  Wj(t  -  kT1  +  r3)  w2(t  -  jT2  +  r4)^ 


Rl212(rl'r2'r3’r4>  =  R1 122  ^  rl ’ tVt2'ta) 

as  the  summation  can  be  rearranged  to  have  an  identical  appearance  to  Rjj22; 

but  with  the  ordering  of  the  r^'s  changed.  Conjugation  can  be  adjusted  to 

convenience  as  the  chip  functions  are  real.  Continuing,  the  remaining 
correlation  products  are: 

R 1 22 1  ^  rl  * r2 ’ r3 ’  r4^  =  R1122^rl 1  r4’  r3’  r2^ 

R2 1 12  ^  rl’  r2’  TV  r4^  =  R1 1 22  ^  r3  ’  r2  ’  rl !  r4^ 

R2 12 1  rr  r2  ’  r3'  r4^  =  R 1 1 2 2 ^ r4 ’  T2,r3,T0 

R2211^rl,r2,r3’r4^  =  R1122^r3’  r4?  rl ’  r2^ 


The  interesting  term  in  the  expression  for  ^122^  rl’ r‘>’ r3’ r4^  KS  t^ie 
complex  exponential: 

n  tT  r  T  T  )  -  J-  e“iu/  ^r2  '  rl)  eiu,,,(r4  "  r3) 

*Wr2’r3’  4^  ’  75r  e  1  1  e  4  J  ’ 

as  the  spectrum  is  obtained  from  the  autocorrelation  by  the  transform; 

00 

SU)  =  /  R(r)  e'iu;r  dr, 

and  only  the  complex  exponentials  depend  upon  r.  The  different  terms  develop 
different  b  functions  that  collapse  the  autocorrelation  products  into 
partial  spectrums  that  may  be  easily  interpreted. 

Evaluation  of  the  function  is  most  easily  performed  by  considering 
the  integral; 

S1234(^  =  /  D^rl’r2’r3’r4)  dr 

-x 

X 

=  jL  f  r2  -  rl>  eiu"(  t4  ~  T3>  e~ia,T  dr 

-X 

For  the  various  correlation  products,  the  corresponding  delay  differences 
can  be  tabulated  in  terms  of  the  ordering  of  the  arguments  for  •’ 


1122 

(r,,r2,r3,r4) 

(r2  ~ 

rl)’ 

(r4  - 

r3^ 

rd,  rd 

2211 

(rvr4,rvT2) 

(r4  - 

r3), 

(f2  - 

rl) 

rd’  rd 

1212 

( T 1  ’  r3  ’  r2  ’  r4  > 

<=> 

(r3  ' 

rl)’ 

(r4  ~ 

Tci) 

<=> 

T,  T 

1221 

f  r, ,  r4>  r3'  rj) 

(  r4  ~ 

rl), 

^  r2  ” 

r3> 

rd  +  r’ 

2112 

<=> 

(r3,r2,r,,r4) 

(r2  “ 

r3h 

<r4  - 

rl) 

<=> 

Td  ~  r’ 

■2121 

<=> 

<Vr2’r3’rl) 

(r2  ~ 

T<0  ' 

(rl  “ 

-r,  -r 

where  the  actual  values  of  the  arguments  of  the  correlation  products  have 
been  used  on  the  right  of  the  table. 


With  those  pairs  the  integrals  become: 

c  -i(w'  - 

S1122M  =  S(u)  e 

Sooiil-')  =  ^1122^'^ 

S 12 i2 ('*'’)  =  6^'  "  +  ^ 


S1221^ 


=  £  /  e-iw'(rd  +  eiw^rd  “  r)  e’1 


-1 UT 


=  e-i(^’  -  w")r 


d  a;*  +  u>"  +  uj) 


=  ei(2tJ"  +  ^)rd  6(u'  +  u>"  +  u;) 


S2112^ 


a-iW  (rH  -  t)  gia;"(rd  +  r)  c-iu;r 


e  1U"  dr 


=  e-1^'  "  w")rd  6(u'  +  w"  -  w) 
=  ei(2w"  "  Td  $(w'  +  a;"  -  u) 


SJ121M  =  *(w'  +  J*  +  W) 

Lastly,  the  autocorrelation  product  spectra  can  be  calculated  from 
these  results.  The  delta  functions  are  inserted  and  the  integrals 
contracted,  when  possible.  The  partial  spectra  are: 


DC  Terms: 


Si  190(0;) 


2 it  6(u)  A2A2  do/'  |V1(u>')|2 

1  -00 


-no  t  , 
e  d 


n  rj  <w’  lv2<“")l2 


iw"r. 


Now  from  Appendix  II,  each  of  the  integrals  is  simply  the  convolution  of  the 
chip  shape  with  itself  evaluated  at  =  Tc/2.  As  this  is  a  unit  height 


triangle  with  a  base  of  2T..  the  result.  is  simply  1/2  if  the  ehip  durnt 
are  approximately  equal. 

Therefore, 


S1122^^  =  2n  ^  ?  A1A2 
=  S22n(u;). 


Chip  Mixing  Terms: 


00 

Sl2l2(o;)  =  k\k\  ^  y~  f  lwi(w"  - 

1  z  -co 

which  will  subsequently  be  evaluated. 


w)I2|v*k)|2. 


s2mM  -  A2Aj  jj  r  r  /  dw"  |W,(-  u”  - 

1  L  -00 


U)|2  |V*(u.")| 


=  S1919(u;) , 


upon  letting  u"  =  a;"  and  using  the  fact  that  the  chip  shapes  are  real. 


c  f  \  a1UTi  .2,2  111  f 

S1221(cj)  =  e  A1A2  2i  T7  TT  J  ' 

1  Z  -oo 


J  dw"  |Vj(w"  *  a>)|2  |vj(w")| 


2  i2a,"rd 
e 


S2112(“>  =  e  1U,r<l  A1A2  nfr/  du"  |Wl(w"  *  w)|2  lW2(“">l2  e 

1  2  -00 

* 

=  si991  (u>) , 


-i2JV. 


md  the  sum  S^i^)  +  ^2112^  's  ro;i*  • 


Evaluation  of  the  Integral: 


All  four  terms  involve  the  integral: 


X 

I(u>,r)  =  ’p’  J*  dtJ"  +  a’)l^  |W2(-”)|2 


-i2u/'r 


1  2  • 


where  r  is  (0,r^).  To  evaluate  the  integral  define: 


F  (u>)  £  /  |V*(u/')|2,  p  =  1,2. 
v  P  v 

Then  F  (u>)  is  the  Fourier  transform  of  the  convolution 
P 


f00  * 

fD(t)  =  J  du  wp(u)wp(u  +  r) 


P  -00 


which  is  the  usual  triangle  of  unit  height  for  rectangular  chips  wp(u) 


Rewriting  the  integral  in  terms  of  the  Fp's: 


f®  * 

I(u>,r)  =  J  du/'  F^u/'  +  u)  F2(u;")  e 


-i2u/'r 


«  r00  r”  f00  iu/'(u  -  v  -  2r)  ion  * 

=  T5r  J  da;"J  duJ  dv  e  e  ^i(u) 


“X  -05  “X 


Using  the  delta  function  obtained  by  the  integration  over  a/': 


f°°  * 

I(u>, r)  =  J  du  f*(u)  f 2 ( u  -  2r)  eiu7U, 


which  is  comparatively  simple  to  evaluate  in  the  cases  of  interest. 


^r^-rr-K  T- 


1-k  rrva«v'\  •%.'  ’  ^rr^r*  vuirrr^i  *vwv»'J  w  ttw*tw  v  w\j  wjbwvif  '  ^/vivrwmirjnyiM  r1*  ^fwT^^wwv^u  «  Trwtnmf  r-W'vrvT'-vr^n-^rrv^  - 


r  =  0: 


I(w'.O)  =  j  du  f | ( u)  f2(u)  eiu;u, 
—00 

As  the  chips  are  approximately  equal: 


f i (u)  f 2 ( u )  = 


1 2 


<T,, 

0,  otherwise. 


and 


r°°  * 

l(w,0)  =  J  du  f^(u)  f 
-00 


i(u)  el£i/u, 


l(w,0)  =  T1  J  du  (1  -  u)2  e  TlU  +  f  du  (1  -  u)2  e  ^ 

L  f)  »  n 


-iulLul 


I(.,0)  =  T,  - [l  - 

1  (JTj? 


which  is  always  positive  and  has  a  value  of  4Tj/6  at  DC.  With  this  result 
the  partial  spectrums  Sj2j2(u;)  ^2121^  are  S^ven  by: 


Sl212(o;)  =  4T1  A2A2  — i-j 

(u/Tj) 


sin(o/T1) 

1  jq — 


=  S2121(w). 


I(w’2rd)  =  /  du  f*(u)  f2(u  -  2 T,  )  eluM, 

“00  u 

i ur,  r *  * 

1  ( ^ ’ 2 rd )  =  e  J  du  Vu  +  T-d)  ^2(u  "  rd)  elUA1’ 
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As  the  chips  are  approximately  equal,  and  only  overlap  in  the  region 


-rd<  u  <  Td: 


fl<u  *  rd>  f2<“  “  rd>  = 


1  "  u" 

J  -  ITJ  .  M  <  rd> 

0,  otherwise. 


I(w,2  r.)  =  e 


i„r.  fV2 

d  du 

-T./2 


'i  r  u 


i 


upon  changing  variables  of  integration.  Evaluating  the  integral  gives: 

i u/r,  1  0  fsin(u;T1/2) 

l(u,2rd)  =  Tj/2  e  [^Tj7?]  [ — uTJ1! - cos(u;T1/2) 


and,  therefore, 


S91  io(w)  -  T  M 


T1  ,2.2  r  1  1 2  fsin(o;T1/2) 


T  A1A2  [uTjJl] 


-  cos(o/T1/2) 


Combination  of  the  Chin  Mixing  Terms: 


With  the  results  in  hand  for  the  four  chip  mixing  terms  (I,  II,  III, 
and  IV),  the  addition  of  them  produces  the  total  mixing  spectrum.  The 
partial  spectrums  an(^  ^1221  ^  can  a(^e(*  t0  Pr°duce: 

§91  io(w)  +  $1991  (uj) 


-  89^9(0;)  +  S91 .9(0;) 


99  1  2  sin(  u/T •.  /  2 ) 

=  4T1  A1A2  fcTj]  — ljTTT2 - cos(wT1/2) 


and  the  other  pair  of  mixing  terms  give: 


^1  91  9  M  +  s 

9 1 9 1  ( w) 


=  8T,  khl  — — , 
1  1  2  U’T ,)* 


s  i  n  ( u/T  L ) 


S3 


Finally,  the  chip  mixing  spectrum  is  the  relatively  simple  result: 


M  A 

$12(^0  =  ^2 1 1 2  ( u;) +  ^1221  ^1212^^ +  ^2121^^ 


-  4T  -  — 

41 1  a1a2  ,  T  , 
(o,'Tt ) 


^  [2  -  cos  ( /2 )  ]  1  + 


sin(wT1/2) 


At  the  chip  rate  harmonics  (u  =  =  2 7r/T^ ) ,  where  a  narrowband  output 

filter  can  be  placed  to  detect  the  chip  rate  line,  the  continuous  chip 
mixing  spectrum  has  the  value: 


1  12  ”2 


which  is  a  major  result! 


Sroiwvxa 


SOESCii 


APPENDIX  V 


Chip/Noise  Mixing;  Terms: 


DC  Term 


Mixing  of  noise  with  a  chip  stream  is  typified  by  the  term: 

Rxx.m<r>  *  XX  Etckcj]  "l*1  -  jTl>  V1  -  kTl  +  rd> 

*k  j  \ 

■E[n  (t  +  r)n(t  +  +  r)]y, 

* (A\«(rd)X  ui(t  ■  kTij  ui(t  -  kTi  *  rd>)’ 

k 

which  is  simply  a  DC  term  as  it  is  independent  of  r.  Stationarity  of  the 
noise  process  and  the  second  moment  of  the  chip  stream  were  used  to  arrive 
at  this  result. 

The  chip  summation  is  no  more  than  the  continuous  spectrum  of  the  input 
chip  stream  evaluated  at  (See  Appendix  II) .  For  rd  of  one-half  of  a  chip 

time,  that  autocorrelation  is  .5.  Therefore,  the  DC  noise/chip 
contributions  are: 

RUrm<r>  =  Z  Vrd>  ' 

=  R22nn(rd)’ 

Rnnll^rd^ ’ 

=  Rnn22(7V’ 

As  the  output  contains  all  four  terms,  the  total  DC  contribution  from  the 
chip/noise  mixing  is: 


=  2A2 


J  Rnn<rd> 
—00 


-U1T 


dr 


or 


=  2TRnn(rd)  S(0). 
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Chin /No  iso  Continuous  Snort ral  Tot  ms: 


Evaluation  of  these  terms  proceeds  much  like  earlier  work.  Calculate  a 
single  term  in  general  form  and  the  other  terms  are  special  cases.  Now, 
consider 


Rlnln  =  XX  E[ckCjJ  wl(t  ~  jTl  +  rl)  wl^  "  kTl  +  r3^ 
*k  J 

•E[n  (t  +  r2)n(t  +  r4)]V 


Rlnli.(r)  =  R  nn^r4  "  r2>  Z  “l(t  '  kTl  *  rl)  “!(t  '  kTl  +  r3>)’ 


=  R 


lnln^rl’  TV  r3’ 


lnnl^l’  TV 

r3’ 

r4> 

=  Rlnln^  71 ’ 

r2’ 

r4’  r3> 

nlln^l’  TV 

r3’ 

r4^ 

=  Rlnln^r2’ 

7r 

r3’  r4> 

nlnl  ^  71  ’  TV 

r3’ 

r4^ 

=  Rlnln^r2’ 

Tv 

rr  r3) 

From  earlier  work,  the  periodic  chip  stream  time  average  is  given  by: 

(a2  X  wl^  “  kTl  +  wl^  ”  kTl  +  r3^) 


iu;,  kT< 


■  <*2  Zh  *  1  ‘> 


=  a r 


with 


/•°° 

2  i  r  * 

=  A  J  w^(t)wj(t  +  r3  -  Ti)  dt, 
1  “00 


=  Rvw<r3  -  rl> 


Therefore,  the  general  term  can  be  written  as 


and  as 


RInnl ( rI •  r2’  r3’  rd  =  Rlnln(rl’  TV  V  r3* 


'  Rn„<r3  -  r2>  KJr4  "  rl>' 


Rnlln*rl’  t2’  r3’  r4^  ‘  Rlnlidr2>  rl’  r3’  r4> 


=  R»n(r4  *  rl>  Rw«<r3  -  rl>’ 


Rnlnl<rr  r2’  r3-  Ti>  =  Rlnln(r2-  rl  •  V  r3> 


*  Rnidr3  -  rl>  Rm<r4  '  r2> 


With  these  general  forms  the  actual  values  of  the  r^'s  may  be  inserted  from 
the  table: 


Rl„lr,(rl'  r2’  r3-  r4>  =  Rn„<r>  R«»(r)’ 


ainndrl'  r2-  r3'  r4>  1  R»n<r  '  rd>  Rw»<r  +  Ti>  • 


Rntln(rl-  r2'  r3'  r4>  =  Vr  *  rd>  R-(r  -  rd)> 


ww 


Rnlnl<rl-  r2’  r3>  *  Rnn(r>  R»»<r>' 


Conjugation  of  the  real  noise  autocorrelation  function  has  been  performed  to 
simplify  subsequent  spectrum  calculations. 


The  spectrum  of  the  term  is  given  by  the  transform  on  r: 


W">  ■  /  O  -  rd>  R™<r  *  rd>  ^ 


The  spectrum  is  easiest  to  evaluate  in  this  form  if  assumptions  are  made 
about  the  about  the  autocorrelation  (or  spectrum)  of  the  noise.  If  the 
noise  spectrum  is  assumed  to  be  a  bell-shaped  one,  with  a  bandwidth  equal  to 
the  chip  rate,  the  the  autocorrelation  of  the  noise  may  be  taken  to  be 


Rnn(r)  =  NIP  u(r)’ 


where  u(r)  is  a  triangle  of  a  base  width  of  2Tj.  If  the  noise  is  taken  to 


have  a  wider  bandwidth,  then  the  triangle  is  more  narrow.  With  this  choice 
of  the  autocorrelation: 


rww: 


'  "V  t  'V  VI 


a 


asa 


a 


Vi« 

1«:C 


lCVZVV 


uu 

sinm^’>  =  nip  a2  /  Vr  -  rd>  u(r  *  rd>  e‘ur  dr' 


which  is  identical  to  integrals  evaluated  previously  in  terms  of  chip  mixing 
( I  ( — ■ .  r ) ) .  The  result  for  the  sum  of  the  four  terms  may  be  written 
immediately  from  the  chip  mixing  results. 


4^)  -  Slnln<^  +  Snl„l<w)+  Slnn2^')  *  Snlln<w) 


0  1  [  r  1  sin  (  /2) 

=  4T1  A1NIP  ~T  71  t2  “  cosC^i/2)]  1  +  z  — u-T  72  ” 
(^TJ  L  L  r 


At  the  chip  rate  harmonics  (u>  =  =  2 tt/T^ ) ,  where  a  narrowband  output 

filter  can  be  placed  to  detect  the  chip  rate  line,  the  continuous  chip/noise 
mixing  spectrum  has  the  value: 


S?nH>  ■  T1  AM  7 


which  is  a  major  result! 


Constant  Input  Spectrum  to  the  Chip  Rate 


A  similar  calculation  for  the  rectangular  spectrum  results  in  a  very 
small  value  for  the  value  of  the  chip/noise  spectrum  at  the  chip  rate.  By 
using  the  constant  valued  spectrum  for  the  noise  (to  the  chip  rate)  the 
value  of  the  continuous  chip/noise  spectrum  in  this  instance  is: 


S^CH>  =  T,  k\  Nj  i  [2  Si(4x)  -  Si(8r)]  , 


where  the  function  Si(x)  is  given  by 


ii(x)  =  f 


To  a  good  degree  of  approximation,  the  values  of  Si (47t)  and  Si (8tt)  can  be 
replace  by  2/x,  the  asymptotic  value  of  the  function. 


asymptotic  value  of  the  function. 


Therefore,  the  continuous  spectrum  has  the  value: 


sc/c 

5ln 


■  T!  NI 


1 

V 


for  the  rectangular  input  noise  spectrum.  This  value  is  considerably 
smaller  than  that  for  the  chip  shaped  input  noise  spectrum  result,  above,  as 
the  chip/noise  products  are  less  correlated  due  to  the  wider  noise 
bandwidth. 
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S: 
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APPENDIX  VT 

Analysis  of  the  Chin  Rate  Detector  Output  S ignal-to— Noise  Ratio 


The  signal-to-noise  ratio  for  a  particular  chip  stream  output  with  a 
pair  of  chip  stream  inputs  and  noise  is  (Section  6.1.4): 


SJ  +  S “  +  a  7r  / 4  NJ  +  3[S1S2  +  05^  +  0 S^J 

with  a  and  &  constants  depending  upon  the  noise  spectrum  model  chosen.  The 
other  chip  stream  has  the  same  expression  for  the  signal-to-noise  ratio;  but 
with  the  numerator  identifying  that  chip  amplitude  (denominators  are 
identical) . 

With  a  rectangular  filter  (with  an  upper  cut-off  of  the  chip  rate) 
9  9 

a  =1  and  0  =  ;r  / 6.  For  a  chip  shaped  impulse  response  filter,  a  jr  /4  =  3/2 

and  J  =  1 . 

For  interpretation  of  the  result,  two  cases  are  of  interest;  large  and 
small  input  signal-to-noise  ratios. 

Small  Input  Sienal-to-iNoi 


Given  the  input  signal-to-noise  ratio 

’SI  4  Sl 

ffji,i  "  nj* 

«  i, 

rewrite  the  expression  above  as: 


So  +  a  7r  /4  N 


0S^l  + 


S^/SJ  +  a  X 74  NJ/SJ  +  3[S2/S1  +  /JSgNj/SJ] 

and  the  result  depends  upon  the  ratio  of  the  two  chip  stream  powers 
Case  I:  The  Chin  Stream  S,  Dominates: 


Si  »  S9 


ft 


■llltl 


rs 

l\l 


TV 


o»l  2  \'.  2t  /  p  h 

a  j  /4  Nj/Sj 


1  4 

*7 


TV 


2 

i,r 


which  is  the  well  known  result  for  a  chip  rate  detector  with  only  a  single 
chip  stream  and  noise  in  the  input  (a  =  1). 

Case  II:  The  Chip  Stream  SL  Dominates: 

S2  >>  S1 


and 


rs 


TV 


0,1  s|/s‘^  +  a  tt2/4  N^/sJ  +  3/S2Nj/S 


1 


Now,  things  are  a  good  deal  more  complicated.  However,  the  results  can  also 
be  written  in  terms  of  the  signal-to-noise  ratio  of  the  second  chip  stream. 

Letting  the  second  chip  stream  signal-to-noise  ratio  be: 


S' 

SJi,2 


A  S2 


then . 


[S' 


o,l  = 


1  +  a 


l 

a 

rJT? 

T 


2  rq-l  2 

N Ji,l 


TW 


i  2  +  3  | 

i,2  a 


'2 

2 

rS] 

5 

7T 

i,2 

and  two  subcases  present  themselves  according  to  whether  the  signal-to-noise 
of  the  second  chip  stream  is  small  or  large.  If  the  second  chip  stream  has 
a  small  signal-to-noise  then  the  output  signal-to-noise  ratio  reduces  to 
simply  the  single  chip  stream  case  as  expected. 

For  a  large  s  gnal-to-noise  of  the  second  chip  stream,  the  output 
signal-to-noise  ratk  becomes: 


o,l 


rs, 


n2 


which  is  worse  than  simply  noise  by  the  factor  of  a  t  /4 
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For  a  large  input  signal-to-noise  ratio 


A  S1 


*]i,l  = 


»  1, 


rewrite  the  expression  above  as: 


S' 

N  o,l 


1  +  S22/S'\  +  3S2/S1  +  3/?  S2/S1  [N/S]ia 


which  depends  upon  the  ratio  of  the  power  in  the  two  chip  streams, 
two  cases  present  themselves. 


Again, 


Case  I:  The  Chip  Stream  S.  Dominates: 


In  this  instance  the  apparent  output  is  that  of  the  reference  chip 
stream  by  itself.  Only  the  self-noise  of  the  reference  chip  stream  is 
present.  And, 


S1  >:>  S2 


*o,l 


Case  II:  The  Chip  Stream  So  Dominates: 


S2  S1 


*  o,l  “  71 


The  output  signal-to-noise  ratio  depends  only  upon  the  ratio  of  the  chip 
stream  powers,  as  might  be  expected. 


Calculation  of  the  signal  and  noise  power  for  N  chip  streams  in  the 
chip  rate  detector  input  is  a  relatively  simple  generalization  of  the  two 
chip  stream  case.  From  Section  6,  an  examination  of  the  terms  in  the  chip 
detector  output  reveals  that  there  are  only  two  types  that  contribute  at  the 
chip  rate.  One  is  the  self-noise  of  the  chip  stream/noise,  which  is 
quadratic  in  the  chip/noise  powers.  There  are  N  chip  terms  and  one  noise 
term.  The  other  type  is  the  intermodulation  noise  of  the  chip  streams/noise 
which  results  in  the  cross-power  terms.  There  are  N(N-l)/2  chip/chip  terms 
and  N  chip/noise  terms. 

Therefore,  the  total  "noise"  out  of  the  chip  rate  detector  is: 


TN 


n  =  X sj * a ,2/4 "MIX sjsk *  3/? X sjNr 


J 


J  k 

j*k 


J 


where  the  sums  are  over  the  N  chip  streams.  By  making  reasonable 
assumptions  about  the  distribution  of  the  powers  (signal-to-noise  ratios)  of 
the  chip  stream  inputs  an  average  output  noise  can  be  calculated.  This 
average  over  the  powers  is: 


E [TN]  =  Nj  [n  E[SNR2]  +  a  *2/4  +  |  N(N-l)  E [SNR] 2  +  3 0  N  E [SNR] 


The  output  chip  line  power  (normalized  to  the  input  noise  power)  for  the  j 
chip  stream  is 


-  th 


S  =  S2/N^  TV 


with  the  average  output  equal  to: 


2n 


E[S]  =  E  [SNR  ] 


With  these  quantities,  the  average  output  signal-to-noise  ratio,  over 
all  chip  streams  is: 


SNR. 


E  [SNRZ]  TV 


av 


N  E[SNR2]  +  a  x2/4  +  \  N(N-l)  E [SNR] 2  +  30N  E[SNR] 


1 

N  7 


TV 


1  +  a  ?r/(4N)  +  i  (N-l)  E[SNR]7E[SNRZ]  +  3/?  E[SNR]/E[SNRZ] 


'  ^  '.V 


Further  evaluation  of  the  average  signal-to-noise  ratio  requires 
specification  of  the  distribution  function  of  the  chip  stream 


signal-to-noise  ratios. 


Two, cases  will  be  considered.  One  will  be  a  uniform  distribution  over 
<iB  of  the  signal-to-noise  ratios  of  the  chip  stream  transmitters.  The  other 
will  be  a  uniform  distribution  over  transmitter  ranges,  with  all  transmitted 
powers  equal. 


Case  I: 


Cn i fo rm  Distribution  of  Signal-to-Noise  Ratios  (in  dB) . 


Assume  the  distribution  is  given  between  zero  and  dB  ,  then 


1  fl,  o  <  dB  <  dBmax , 
Pd“(  ’  =  Jo,  otherwise. 


and  obtaining  the  distribution  of  S  through  the  relation 
dB  =  10  log10(S): 


PS(S)  = 


1  <  s  <  S. 

-  max 


The  expectations  follow: 


E[S]  = 

max 


—  i 

1  Jmax  1 


and  the  desired  ratios  are  (for  S  „  >>  1) 

v  max  ' 


EfSl2  =  _ 2 

E[S2]  TF^ 


E  f  SI  = 

E[S2] 


With  these  relations,  the  average  chip  rate  detector  output  signal-to-noise 
ratio  is: 


av  -  N 


1  ,  a  //(4N)  ♦  3(N-l)/ln  *  6/?/SBax 


For  a  large  number  of  emitters  (N  >>  1)  and  high  input  maximum 

signal-to-noise  ratios  (S„v  >>  100  (20  d  B)),  the  the  formula  may  be 

InctX 

approximated  by: 


TV 


1  +  3 ( N— 1 ) /In  S 
v  ' '  max 


A  conclusion  can  be  reached.  The  average  signal-to-noise  ratio 
decreases  at  least  as  fast  as  N  for  emitters  spread  over  extremely  wide 
dynamic  ranges.  For  a  60  dB  dynamic  range: 


SNR 


.  1 
av  N 


TV 


1  +  .22  (N— 1 ) 


Case  II: 

Uniform  Distribution  of  Ranges 

Assume  the  distribution  is  given  between  R.„  and  R  „  .  then 

°  min  max 


pdB(dB)  =  R-  R  • 
max  mm 


1?  Rmin  -  R  -  Rmax’ 


[0,  otherwise, 
and  obtain  the  distribution  of  S  through  the  relation 


As  in  case  I,  for  large  max  to  min  range  range  ratios: 

TV 


SNR 


1 


av  "  N  7 


1  +  a  //(4N)  +  %  (N-l)  Rm;n/Rmav  +  9/3/SNRm-  rantTO 
1  K  '  2  v  ’  min'  max  '  min  range 


Again,  assuming  the  minimum  range  signal-to-noise  ratio  is  large,  and  N 
also: 


SNR 


av 


1 

N 


TV 


1  + 


J  (N-»  R»iAax 


A  similar  conclusion  to  the  uniform  distribution  in  dB  case  can  be 
reached.  The  average  signal-to-noise  ratio  decreases  as  the  number  of 


VT  V  WJ1 


emitters  (N).  A  choice  of  maximum-to-min imnm  range  of  1000  (60  HR 

dispersion)  gives,  in  this  instance: 


av  a  r  l 

1  +  .Gu45(N-l)j 

Summarizing,  a  large  number  of  emitters  tends  to  obscure  each  emitter's 
chip  rate  line.  The  obscuration  increases  as  10  log(N)  in  dB! 
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NOMOGRAM  DETERMINES  PROBABILITY  OF 


Detecting  Signals  in  Noise 


By  DONALD  E.  BAILEY 
NEIL  C.  RANDALL 

Hfftarch  Division.  Philco  Corporotloo.  Phllfdtlphlt.  Po. 


/  htir  through  points  of  tn>n  l.nuu-u  vanublrt  will  girr  thr  third 
i-rt  i  to  tilr .  wit  h  or  r  n  i  n  c  y  *  nnf  i!  1 1  i  s  in  pi  n  l  t  hr  y  /  pV  im/r.  If  n  S/K  rat 
ip/  s  etb  if  uvmloblr  uuri  o  /'/.  fi/  /O'*  ts  uautrd,  a  straight  Itur  bttu'tem 
thr$i  piipnff  giurm  thr  |/i  ubnbilil  y  nf  rirtrctinn  at  I3J  pnrrnt.  Thrtsh* 
••Id  Itvtl  it  about  i.f 


In  rkceivinc  systems  it  is  often 
desirable  to  detei-mine  detection 
probability  for  signals  in  noise. 
This  nomogram  gives  a  method 
of  doing  it  without  solving  in¬ 
tegral  equations'  ‘.  or  using  ta¬ 
bles*'*.  when  dealing  with  envel¬ 
ope  detection  of  signals  with 
nonfluctuating  amplitude  accom¬ 
panied  by  additive  Gaussian 
noise. 

Signal-to-noise  ratio  S/S  is 
the  rms  power  signal-to-noise 
ratio,  threshold  level  B  ia  the 
input  voltage  above  which  the 
detector  gives  an  output,  and 
false-alarm  probability  P„  is  the 
probability  that  noise  voltage 
will  exceed  the  threshold.  Equ*^ 
tion  1  relates  threshold  level  and 
false-alarm  probability,  assum¬ 
ing  Gaussian  noise  and  normal¬ 
izing  with  respect  to  rms  noise 
current. 

Probability  P,  that  the  sum  of 
instantaneous  voltages  of  signal 
and  noise  will  exceed  the  ti.r«M*.- 
old.  is  the  integral,  from  thresh- 
old  level  to  infinity.*' of  the  prob- 
ability  density  function  of  signal 
plus  noise1  \  as  in  Eq.  2.  where 
J.  is  the  Bessel  function  with 
imaginary  argument,  and  V  the 
envelope  of  the  sigi'.al-pliis-noise 
current  normalized  by  rn.*  iw.ivi- 
current.  Density  functions  m  the 
figure  show  typical  rehit  iuo'h  i  |i> 
among  P,..  B  and  i>*/.V.  The 

nomogram  is  n  semivmp:  rical 
solution  of  Eq.  2  for  a  pwwer 
S/S  ratio  measured  in  the  re¬ 
ceiver  i-f. 
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